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Astrophysical black hole candidates are thought to be the Kerr black holes of general 
relativity, but there is not yet direct observational evidence that the spacetime geometry 
around these objects is described by the Kerr solution. The study of the properties of the 
electromagnetic radiation emitted by gas or stars orbiting these objects can potentially 
test the Kerr black hole hypothesis. This paper reviews the state of the art of this 
research field, describing the possible approaches to test the Kerr metric with current 
and future observational facilities and discussing current constraints. 
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A. Motivations 


Einstein’s theory of general relativity was proposed in 


1915 (Einstein, 1916). The first experimental confirma¬ 


tion was attained in 1919, from the measurement of light 
bending in the vicinity of the surface of the Sun by Ed¬ 


dington (Dyson et al. 1920). Thanks to this observation, 


the theory immediately became very popular. However, 
systematic tests of general relativity started much later. 
Experiments in the Solar System started in the 1960s. 
Tests using the observations of binary pulsars began in 
the 1970s. In the past 50 years, a large number of exper¬ 
iments have confirmed the predictions of general relativ¬ 


ity in weak gravitational fields (Will 2014). The focus of 


current experiments has now shifted to testing the theory 
in other regimes. 

Tests of general relativity at large scales are mainly 
motivated by cosmological observations. While the dark 
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matter problem is more likely due to some weakly in¬ 
teracting massive particles beyond the Standard Model 
of particle physics rather than a breakdown of Newton’s 


law of gravitation (Bergstrom 20001, see in particular 


the Bullet Cluster (Clowe et al. 


2004), the problem of 


dark energy is completely open (Frieman et al. 2008). 
At the moment, the accelerating expansion rate of the 
Universe may be explained by a small positive ad hoc 
cosmological constant. However, it is possible that the 
actual explanation is either a breakdown of general rel¬ 
ativity at large scales or the existence of some new field 
with peculiar properties. In the past 20 years, testing 
general relativity at cosmological scales has been a very 


active research field (Daniel et al. 2010). 


More recently, there has been increasing interest in 
testing general relativity in strong gravitational fields, 
which is the regime where the theory is more likely to 
encounter deviations from its prediction^ One of the 
most fascinating predictions of general relativity is the 
existence of Kerr black holes (BHs). The so-called BH 
candidates are astrophysical compact objects that can be 
naturally interpreted as the Kerr BHs of general relativ¬ 
ity or they could be something else in the light of new 
physics. 

A direct observational confirmation of the nature of 
these objects could be seen as an important test of 


a reliable and sophisticated astrophysical model to test 
fundamental physics. Note that such a research field is 
young and still a developing one. While tests of general 
relativity in the Solar System and with binary pulsars 
started about 60 years ago, studies to directly test the 
nature of BH candidates are recent and became a hot 
topic only in the past 10 years. As a not yet fully mature 
research field, in order to make progress it is common to 
employ some approaches lacking the necessary scientific 
rigor and that are sometimes criticized. 

The results of current electromagnetic observations can 
be summarized as follows. BH candidates cannot be as¬ 
trophysical bodies with a surface made of ordinary mat¬ 
ter (jBrodericket^ 2009 Narayan & McClintock , 2008), 
nor can be compact stars made of exotic non-interacting 


particles (Bambi & Malafarina 2013). Even if the con¬ 


straints are currently weak, the quadrupole moment of 
these objects seems to be in agreement at the level of 


30% with that which is expected for a Kerr BH (Val- 


tonen et al. 2010). The thermal spectrum of the ac¬ 


cretion disk around the stellar-mass BH candidates in 
GRS 1915-1-105 and Cygnus X-1 looks like that expected 
for a very fast-rotating Kerr BH, and this can constrain 
at least some kinds of deviations from the Kerr geom¬ 
etry (Bambi 2014a Kong et al. 2014). Similar con¬ 


general relativity in the strong gravity regime (Bambi 


straints can be obtained for supermassive BH candidates 
under some reasonable assumptions about the evolution 


2011a 2013a). Deviations from the Kerr metric may of the spin parameter (Bambi 2011b). These results are 


be expected from classical extensions of general rela¬ 


tivity (Berti et al. 2015 Capozziello & de Laurentis 


2011) as well as from macroscopic quantum gravity ef¬ 
fects ( jDvali fc Gomez 2013a|b Giddings 2014). Indeed, 


there are some arguments suggesting that the critical 
scale at which classical predictions may break down is 
not at the Planck one, as expected in a collision between 
two particles, but at the gravitational radius of the sys¬ 
tem (Dvali & Gomez, 2013a|b Giddings, 2014). 


reviewed in the next sections. However, there is much 
work to do to test the nature of BH candidates, and new 
and more powerful observational facilities are necessary 
if we want to achieve this goal. 

Gravitational wave tests have been reviewed, e.g., 
in Yunes & Siemens (2013). Gravitational wave con¬ 


There are two possible approaches to test BH candi¬ 
dates: electromagnetic and gravitational wave observa¬ 
tions. The aim of this paper is to review current attempts 


straints on alternative theories of gravity can be obtained 
from the observed decay of the orbital period in binary 
pulsars, see e.g. Yunes & Hughes (2010|. The best system 


for gravitational wave tests may be an extreme mass ra¬ 
tio inspiral (commonly called EMRI), in which a stellar- 
mass compact object slowly falls onto a supermassive 


with electromagnetic radiation. Generally speaking, this 

BH candidate (Barack & Gutler 2007| 

Glampedakis & 

goal is very challenging, because it is difficult to have 

Babak 

2006 

Ryan 

1995 

). Accurate measurements of 


BH quasi-normal modes may also do the job (Berti et 


al. 2006 Gossan et al. 2012). The detection of gravita- 


^ There is no single definition of “strong field” in general relativ¬ 
ity. Even binary pulsar tests are sometimes classified as tests 
in strong gravitational fields. If we consider the strength of the 
Newtonian potential (when the metric is close to Minkowski), 
binary pulsar measurements test weak gravitational fields, be¬ 
cause M/r ~ 10“®, where M is the mass of the system and 
r is the distance between the pulsar and the companion. In 
this context, strong gravitational fields have M/r ~ 0.1-1. A 
more coordinate-independent measure is given by some curvature 
invariant, like the Kretschmann scalar However, 

this introduces a length scale. If we identify the Planck length 
Lpi = ~ 10“®^ cm as the natural scale in general rela¬ 

tivity, no astronomical observation would probably ever reach 
strong fields. 


tional waves by LIGO in September 2015 has opened a 
new window (Abbott et al., 2016a). Gravitational wave 


detectors now promise precision tests of the strong grav¬ 
ity regime in 5-10 years. 

The content of the paper is as follows: an introduction 
section devoted to briefly review the motivations to test 
the Kerr metric, a discussion of some important prop¬ 
erties of the Kerr metric, and a discussion of basic as¬ 
tronomical observations of astrophysical BH candidates. 
In Section |TTj I discuss the general approach to test the 
Kerr metric with electromagnetic radiation, its limita¬ 
tions, and some important phenomena that may show up 
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in non-Kerr metrics. In Section III, I review the possi¬ 
bility of testing the Kerr metric with X-ray observations. 
This is currently the only available electromagnetic ap¬ 
proach to probe the strong gravity field close to these 
objects. Section |IV] is reserved to discuss possible mea¬ 
surements of SgrA*, the supermassive BH candidate at 
the center of our Galaxy, which can be considered as a 
special case. In Section |V) I briefly introduce a few more 
approaches to test the Kerr metric. Section |Vl| discusses 
the differences between tests with electromagnetic and 
gravitational wave observations. Summary and conclu¬ 
sions are reported in Section [Vll| 

Throughout the paper, I employ the convention of a 
metric with signature (—1-++) and units in which Gn = 
c = I, unless stated otherwise. In these units, the size of 
a stellar-mass BH is 


" = (u^) 

In the case of supermassive BHs with mass M = 10® Mq 
(10® Mq), we have M « 1.5 • 10® km (1.5 • 10® km) and 
» 5 s (1.4 hr). 


B. Kerr black holes 


Technically, a BH in an asymptotically flat spacetime 
Ml is defined as the set of events that do not belong to 
the causal past of future null infinity J“(I^), i.e. 


B = M-J-{T+) ^ 


( 2 ) 


The event horizon is the boundary of the BH region. See, 
e.g., Misner et al. (1973), [Wald (1984) or Poisson (2004) 
for the details. In other words, all future-directed curves 
(either time-like or null) starting from the region B fail 
to reach null infinity . A BH is thus an actual one-way 
membrane: if something crosses the event horizon it can 
no longer send any signal to the exterior region. 

While the event horizon is a global property of an en¬ 
tire spacetime, the apparent horizon is a local property 
and is slicing-dependent. The apparent horizon can be 
defined after introducing the concept of a trapped sur¬ 
face. A trapped surface is a 2-dimensional surface S with 
the property that the expansion of ingoing and outgoing 
congruences of null geodesics orthogonal to S is nega¬ 
tive everywhere on S. An apparent horizon is the 2- 
dimensional intersection of the 3-dimensional boundary 
of the region of the spacetime that contains a trapped sur¬ 
face with a space-like hyper-surface. Physically speaking, 
outward-pointing light rays behind an apparent horizon 
actually move inward and therefore they cannot cross 
the apparent horizon. Under certain assumptions, the 
existence of an apparent horizon implies that the slice 
contains an event horizon; the converse may not be 


true (Wald & Iyer, 1991). More details can be found, e.g.. 


Misner et al. (1973), 


Wald (1984) or Poisson (2004). 


Note that there are scenarios beyond classical general 
relativity in which a collapsing object does not gener¬ 
ate a BH with a central spacetime singularity. For an 
observer in the asymptotic flat region, the collapse may 
generate an apparent horizon for a finite time; the lat¬ 
ter may be interpreted as an event horizon if the ob¬ 
servational timescale is shorter than the lifetime of the 


apparent horizon (Bambi et al. 2014a 2016b Frolov & 


Vilkovisky, 1981). 


In general relativity, the simplest BH solution is the 
Schwarzschild metric that describes an uncharged and 
non-rotating BH in vacuum and has only one param¬ 
eter, the BH mass M. If the BH additionally has a 
non-vanishing electric charge Q, we have the Reissner- 
Nordstrom solution. If the BH has a non-vanishing spin 
angular momentum J, we have the Kerr solution. The 
Kerr-Newman metric describes a rotating BH with a non¬ 
vanishing electric charge. 

The fact that BHs are described by a small number of 
parameters (M, J, Q) is explained by the so-called “no- 
hair” theorems ^ wh ich were pioneered in Carter (1971) 
and Robinson (1975) and the final version is still a work in 
progress (Chrusciel et al. 2012). The name no-hair indi¬ 
cates that BHs have no features (hairs), in the sense that 
they are completely specified by a small number of pa¬ 
rameters. In the context of tests of general relativity and 
of the Kerr metric, the “uniqueness” of the Kerr-Newman 
solution is also relevant: under the same assumptions of 
the no-hair theorems, there are only Kerr-Newman BHs. 

“Hairy” BHs generically arise when gravity couples 


to non-Abelian gauge fields (Volkov & Gal’tsov, 1989 


1999) or when scalar fields non-minimally couple to 


gravity, e.g. a dilation field in Einstein-dilaton-Gauss- 
Bonnet gravity (Mignemi & Stewart 1993). An im¬ 
portant example of the violation of the no-hair theo¬ 
rem was presented in (Herdeiro & Radu 2014), in which 
hairy BHs are possible by combining rotation with a har¬ 
monic time-dependence in the scalar field. Another in¬ 
teresting example is Chern-Simons dynamical gravity. In 
this framework, non-rotating BHs are described by the 
Schwarzschild solution, but rotating BHs are different 


from those of Kerr ( 

Yunes & Pretorius 

2009 

). A review 

on hairy BHs is ( 

Herdeiro & Radu 2015 

). A 

review on a 


number of BHs in alternative theories of gravity is (Berti 


et al., 2015). 


If BHs were merely theoretical solutions of the Ein¬ 
stein equations, they would not be so interesting. They 


^ Note that there are a number of assumptions concerning these 
theorems. Specifically, the spacetime must be stationary, asymp¬ 
totically flat, and have 4 dimensions; the only stress-energy ten¬ 
sor is due to the electromagnetic field; the exterior region must 
be regular (no naked singularities) JChrusciel et al.||2012||. 
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become interesting because we have compelling observa¬ 
tional evidence of their existence in the Universe. BHs 
can be created by the gravitational collapse of matter. 
The simplest example for which there is an analytic so¬ 
lution is the collapse of a spherically symmetric cloud of 
dust, the so-called Oppenheimer-Snyder model (Oppen- 


heimer & Snyder 19391. This example shows how the 


final product of collapse is a Schwarzschild BH with a 
central singularity. Scenarios closer to reality require nu¬ 


merical calculations (Baiotti et ah, 2005 Baiotti & Rez 


zolla, 2006). For a general review on gravitational col¬ 
lapse, see e.g. Joshi & Malafarina (2011). A BH can be 


created by the gravitational collapse of a very heavy star, 
after the latter has exhausted all its nuclear fuel and in 
the case the degenerate neutron pressure cannot balance 
the gravitational force. The process should be so com¬ 
mon that from stellar population arguments we expect 
that today there are about 10®-10® BHs in our Galaxy 


formed from core-collapse (Timmes et ah, 1996). 


It is worth stressing that the stationary Kerr solution 
of general relativity should well be capable of describ¬ 
ing the spacetime around astrophysical BHs formed from 
gravitational collapse. In general relativity, initial devi¬ 
ations from the Kerr metric are quickly radiated away 


through the emission of gravitational waves (Price, 19721. 
For macroscopic BHs, the equilibrium electric charge is 


completely irrelevant in the spacetime geometry (Bambi 


et al. 2009). The presence of an accretion disk is nor¬ 


mally negligible, because the disk is extended and has a 
low density (Barausse et al., 2014). Moreover, the disk 


mass is many orders of magnitude smaller than that of 
the BH, so its impact on the measurement of the met¬ 


ric would be extremely small in any case (Bambi et al. 


2014b). 


C. Basic properties of the Kerr metric 

In Boyer-Lindquist coordinates, the line element of the 


Kerr metric reads (Chandrasekhar 1985 Kerr, 1963) 


, 2 2Mr\ , 2 ‘iMars\T?9 , ,, S , 2 

= - 1 -— dt^ - - --- 'didcl) + —dr^ 

\ 2-; y 2 j za 

^ in-? ( ? ? 2Ma?rs\T?9\ , 2 n 1 ,? /o\ 

+Sd9^ -I- f -I- -I- — -j sm^ 9d(j) , (3) 

where a = J/M, E = -|- cos^ 9, and A = — 2Mr + 

a^. The dimensionless spin parametei0is a* = a/M = 
J/M^. 


^ Throughout this review article, I use a* to indicate the dimen¬ 
sionless spin parameter of the spacetime, namely a* = J/M^. 
In the Kerr metric and in the non-Kerr metrics discussed in this 
paper, a/M coincide with J/M"^, but this is not universally true 
when we consider deviations from the Kerr solution. J/M"^ is 
a physically measurable quantity. In the general case, a/M is 


In the Kerr metric in Boyer-Lindquist coordinates, the 


event horizon is defined by the larger root of = 0 (Mis- 


ner et al. 1973 Poisson 20041. The solution of g'^'^ = 0, 


which is equivalent to A = 0, is 

= M zb . 


( 4 ) 


r_|_ is the radius of the event horizon and ranges from 
2M, when a* = 0, to M, when |a*| = 1. r_ is the radius 
of the inner horizon that turns out to be a Cauchy hori¬ 


zon. Metric perturbations are unstable for r < r_ (Droz 
1997 Poisson & Israel, 1990|. This means that the inte¬ 


rior solution of the Kerr metric breaks down for r < r_. 
Deviations from the Kerr geometry should thus be ex¬ 
pected already at r « r_, which coincides with r_|_ in the 
case of an extremal Kerr BH with |a*| = 1, and not just 
at a Planck length distance from the central singularity. 
The exterior solution is regular (no spacetime singulari¬ 
ties and closed time-like curves) only for |a*| < 1. When 
I a* I > 1, there is no horizon and the spacetime describes 
a naked singularity, see Subsection |H.C| 

The properties of the orbits around a BH are an impor¬ 
tant tool to connect possible observational effects with 
the geometry of the spacetime ( {Bardeen et al. 1972). 
Time-like circular orbits in the equatorial plane play a 
special role, because accretion disk models usually as¬ 
sume that the disk is in the equatorial plane and parallel 


to the BH spin (see Subsection HI. A and reference therein 
for more details). 

Circular orbits in the equatorial plane exist only for 
radii larger than a critical radius, called the radius of the 


photon orbit (Bardeen et al. 1972) 


r,y = 2M <1-1- cos 


arccos 


(^s) 


( 5 ) 


where here and in what follows the upper sign refers to 
corotating orbits (orbital angular momentum parallel to 
the BH spin), while the lower sign is for counterrotating 
orbits (orbital angular momentum antiparallel to the BH 
spin). No circular orbits exist for r < r^. Massive parti¬ 
cles can be in the photon orbit in the limit in which they 
have infinite energy. 

The marginally bound circular or6zt separates unbound 
circular orbits (r < rmb) from bound circular orbits (r > 
rmb)- The specific energy of a test-particle is A > 1 in 
unbound orbits (the particle has the energy to escape to 
infinity) and A < 1 in bound orbits. In the Kerr metric in 
Boyer-Lindquist coordinates, the radius of the marginally 


bound circular orbit is (Bardeen et al. 1972) 


Cinb = 2M =F a -I- 2 a/m (M a). 


( 6 ) 


a parameter of a metric to be indirectly inferred, and affecting 
more than just the asymptotical mass-current dipole moment of 
the spacetime. 
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FIG. 1 Radial coordinates of the event horizon r+, the photon 
circnlar orbit r^, the marginally bound circular orbit rmb, and 
the ISCO risco in the Kerr metric as functions of the spin 
parameter a,. For every radius, the upper curve refers to the 
counterrotating orbits, and the lower curve to the corotating 
orbits. 


indicate very compact and massive objects, whose prop¬ 
erties match those of BHs in general relativity. A BH 
is currently the least exotic hypothesis for most of these 
objects. However, many others adopt a less conserva¬ 
tive perspective and call them BHs when there is a dy¬ 
namical measurement of their mass (without testing the 
Kerr-ness of the metric or the possible existence of some 
kind of horizon), bestowing the term BH candidate to 
those objects for which there is no dynamical measure¬ 
ment of their mass but share common features observed 
in sources with a BH. 

BH candidates are grouped into two classes: stellar- 
mass and supermassive BH candidates. There is proba¬ 
bly a third class of objects, intermediate-mass BH candi¬ 
dates (Coleman Miller & Colbert, 2004), hlling the mass 
gap between the stellar-mass and the supermassive ones, 
but the measurement of these objects is more difficult. 
They will not be considered in the following discussion, 
because they are currently unsuitable sources to test the 
Kerr metric with electromagnetic radiation. 


1. Stellar-mass black hole candidates 


Finally, we address the topic of stable orbits. The crit¬ 
ical radius separating unstable and stable circular orbits 
is the radius of the marginally stable circular orbit, rms, 
more often called the radius of the innermost stable cir¬ 
cular orbit (ISCO), risco- In the Kerr metric in Boyer- 


Lindquist coordinates, we have (Bardeen et ah, 19721 


nsco = 


3+ Z 2 T V(3-^i)(3 + Zi+2Z2) 

= 1 + (1 - [(1 + a *)^/^ + (1 - a *) 

^2 = ^3al+Zl. 


M, 

1/3 


(7) 


In the Kerr metric, equatorial circular orbits are always 
vertically stable, while they are radially unstable for r < 
nsco- 

Fig. shows the radial coordinates of r+, rmb, and 
I'lSCO as functions of the BH spin parameter a*. For ev¬ 
ery radius, the upper curve refers to the counterrotating 
orbits, and the lower curve to the corotating orbits. More 
details on the properties of the orbits in the Kerr metric 


can be found in Bardeen et al. (19721. As briefly reviewed 


in Subsection H.B in non-Kerr backgrounds there may 
be substantial differences, which may lead to specific ob¬ 
servational signatures. 


In general relativity, the maximum mass for a compact 
star made of neutrons, mesons, or quarks for plausible 


matter equations of state is 

about 3 Mq (Kalogera & 

Baym 

1996 

Lattimer 

2012 

Rhoades & Rufhni 1974||‘^ 


Within standard physics, if an object is compact and ex¬ 
ceeds 3 Mq it can only be a BH and is classihed as a 
BH candidate. In this context, a reliable mass measure¬ 
ment is thus crucial. BH candidates with a mass ranging 
from 3 to about 100 Mq can be called stellar-mass BH 
candidates. 

Almost all the currently known stellar-mass BH can¬ 
didates are in X-ray binaries. Their masses range from 


about 5 to about 20 Mq (Remillard & McClintock, 20061. 
We can infer that they are compact from their short 
timescale variability. Their mass can instead be mea¬ 
sured by studying the orbital motion of the companion 
star, typically with optical observations. These systems 
can be studied within the framework of Newtonian me¬ 
chanics, because the companion star is always far from 
the BH candidat^^ The key-quantity is the mass func¬ 


tion (Casares & Jonker 20141 


27rGN (M -k M^f ’ 

where Kc = Vc sin i is the maximum line-of-sight Doppler 
velocity of the companion star, Vc is the velocity of the 


D. Black hole candidates 


Astronomical observations have discovered a number 
of “BH candidates” (Narayan 20051. In this article, I 
adopt quite a conservative attitude and I use this term to 


^ In alternative theories of gravity, this bound may somewhat 
change l| Astashenok et al.| |2014| |Doneva et al.| |2013|l . 

® Systems m which the companion star is very close to the BH 
candidate may not exist, because the strong gravitational field 
around the compact object would disrupt an ordinary star. 
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GRS 1009-45 GRS 1124-683 


GS 2000+25 H1705-250 

• • •• 

A0620-00 GRO J0422+32 


GRO J1655-40 

m • 



GS 2023+338 GS 1354-64 


XTE J1550-564 

mm 

GX 339-4 4U 1543-47 


FIG. 2 Sketches of 22 binaries with a stellar-mass BH candidate confirmed by dynamical measurements. For every system, the 
BH accretion disk is on the left and the companion star is on the right. The orientation of the disks indicates the inclination 
angles i of the binaries. The distorted shapes of the stellar companions are due to the gravitational fields of the BH candidates. 
The size of the latter should be about 50 km, compared with the distance of the Sun to Mercury of about 50 million km and 
the radius of the Sun of 0.7 million km (top left corner). Figure courtesy of Jerome Orosz. 


companion star, i is the inclination angle of the orbital 
plane with respect to the line of sight of the observer, 
T is the orbital period, and M and are the mass of 
the compact object and of the stellar companion, respec¬ 
tively. Here I reintroduced Gn because this formula is 
usually presented in this form. If we can get an indepen¬ 
dent estimate of i and and we can measure K^. and T, 
it is possible to determine the mass of the compact object 
M. Note that T and Kc, and therefore the measurement 
of the mass function, can be obtained from light curves 
and spectroscopy. Moreover, from the right-hand side in 
Eq. ([^, we can see that f{M) < M; that is, from the es¬ 
timate of the mass function we can directly infer a lower 


bound on the mass of the dark object. When the mass M 
of a compact object exceeds 3 Mq, the object is classified 
as a BH candidates 

Today we have 24 “confirmed” BH candidates in X-ray 
binaries, where the term confirmed is used to indicate 
that there are dynamical measurements of their masses 
and they exceed the 3 Mq bound. Actually these masses 
are typically higher than 5 Mq, and in some cases the 
same mass function f{M) exceeds 5 Mq, which means 


Massive stars exceeding the bound 3 Mq are not compact objects 
and they can be identified by optical observations. 
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the mass of the compact object is larger than 5 Mq in¬ 
dependent of the estimate of the mass of the companion 
star Me and the viewing angle of the orbit i [because 
M > /(M)]. Fig. [^is a schematic diagram of 22 binaries 
with a dynamically confirmed BH candidate. In addition 
to the 24 confirmed BH candidates, we currently know 
another ^40 BH candidates in X-ray binaries that are 
without a dynamical measurement of their masses. They 
are nevertheless classified as BH candidates since they 
have features of BH candidates, but it is possible that 
some of them are actually neutron stars. 

The presently known stellar-mass BH candidates in bi¬ 
naries are mainly in our Galaxy, and a minor number are 
in nearby galaxies. Note that we expect a population of 


about 10®-10® stellar-mass BHs in our Galaxy (Timmes 


et al. 1996). So we presumably know only a very small 


fraction of the stellar-mass BHs in our Galaxy and this 
is because the detection of these objects is only possible 
under some special conditions related to the properties 
of the binary system. A list of confirmed and uncon¬ 
firmed stellar-mass BH candidates in binaries with their 
properties can be found in 


Ozel et al. 


( 2010 ). 


X-ray binaries can be grouped into 2 classes: low-mass 
X-ray binaries and high-mass X-ray binaries. Low-mass 
X-ray binaries are systems in which the stellar companion 
is not more than a few Solar masses 3 Mq) and the 


mass transfer occurs by Roche lobe overflow (Savonije 


1978). These systems are typically “transient” X-ray 


sources because the mass transfer is not continuous. For 
instance, a similar system may be bright for some months 
and then be in a quiescent state for decades. High-mass 
X-ray binaries are systems in which the stellar compan¬ 
ion is massive 10 Mq) and the mass transfer from the 
companion star to the BH is due to the wind from the 
former. These systems are “persistent” X-ray sources. 
Among the 24 confirmed BH candidates in X-ray bi¬ 
naries, 6 are persistent sources (Gyg X-1, LMG X-1, 
LMC X-3, M33 X-7, NGG 300-1, and IG 10 X-1) and 
the other 18 sources are transient. 

In addition to BHs in binary systems with ordinary 
stars, we expect also the existence of isolated BHs and 
BHs in binary systems with neutron stars or other BHs. 
Attempts to identify isolated BHs are based on microlens- 
ing techniques (Agol et al. 2002 Gould 20001, but these 
observations are difficult and there are currently only 
some weak candidates. BHs in binary systems with neu¬ 
tron stars or other BHs can potentially be discovered with 
gravitational waves in the next years. The LIGO/Virgo 
collaboration recently announced the detection of grav¬ 
itational waves from a binary system in which two BHs 
of about 30 Mq merged into a BH of ^ 60 Mq (Ab¬ 


bott et al. 2016aI. The detection of another binary BH 


merger was reported in Abbott et al. (2016b). Gravi¬ 
tational wave detectors now promise to discover tens of 
stellar-mass BH candidates in the next few years. For 
the time being, we do not know any BH candidate in a 


binary in which the stellar companion is a pulsar, but a 
similar system is not expected to be too rare and there 
are efforts to find BH binaries with pulsar companions 
with radio observations. 


2. Supermassive blacl< Inole candidates 


Supermassive BH candidates have a mass in the range 
M ^ 10® — 10^° Mq (Kormendy & Richstone, 1995). 
They are harbored in the center of galaxies, and it is 
supposed that any normal galaxy has a supermassive BH 
at its center, namely any galaxy that is not too smalQ 
Small galaxies usually do not have a supermassive BH. 
Once again, these objects are classified as BH candidates 
because, from the estimate of their mass and volume, a 
BH is the least exotic hypothesis for most of these ob¬ 
jects. In general, their mass can be dynamically mea¬ 
sured by the motion of the orbiting gas. The estimate of 
the volume follows from the timescale variability of these 
sources. 

The strongest evidence for the existence of a supermas¬ 
sive BH comes from the center of our Galaxy. Much work 
in this direction was done by Genzel and his group (|Gen- 
zel et ah] |1997[ |1996[ [Gillessen et al.[ |2009[ [Schodel ~ 
al. 2002). From the study of the Newtonian motion of 


individual stars, we can infer that the mass of the com¬ 
pact object is about 4 • 10® Mq. The upper bound on 
the size of the BH candidate comes from the minimum 
distance approached by one of these stars, which is less 
than 45 AU and corresponds to 600 M for a 4 • 10® Mq 
BH (Ghez et al. 2005). With such a large mass in a rel¬ 


atively small volume, such an object cannot be a cluster 
of non-luminous compact bodies such as neutron stars. 
This is because the lifetime of the cluster due to evapo¬ 
ration and physical collisions would be shorter than the 


age of this system (Maoz 1998). The most natural in¬ 


terpretation is that it is a BH and, according to general 
relativity, the spacetime geometry around an astrophysi- 
cal BH should be well approximated by the Kerr metric. 


3. Black hole horizon 

The key-feature that makes an object a BH is the ex¬ 
istence of an event horizon. Interestingly, there are some 
observations that show that BH candidates have no “nor¬ 
mal” surface but instead have something that may indeed 
be an event horizon. In particular: 

1. Most X-ray binaries are transient sources: there¬ 
fore they spend a long period in a quiescent state 


^ Exceptions might be possible. For example, the galaxy A2261- 
BCG has a very large mass but it might not have any supermas¬ 
sive BH at its center ([Postman et al.j [2012^ . 

























































with a very low mass accretion rate and luminos¬ 
ity. This is true for X-ray binaries containing either 
neutron stars or BH candidates. However, it turns 
out that BH candidates can be extremely underlu- 
minous in comparison to neutron stars in the quies¬ 
cent state. This should indeed be expected, because 
the thermal energy locked in the gas can be com¬ 
pletely lost when the gas crosses the horizon, while 
in the case of the neutron star the gas hits the sur¬ 
face of the compact object and can subsequently re¬ 


lease energy (Narayan & McClintock, 20081. More 


interestingly, during the quiescent period, neutron 
star binaries show a thermal blackbody-like compo¬ 
nent in the X-ray band, which is interpreted as the 
emission from the neutron star surface. No similar 
component is found in binaries with a BH candi¬ 


date (McClintock et al. 2004 1 . The fact that there 


has been no detection of any thermal component 
may be interpreted as evidence for the absence of 
a normal surface and therefore the existence of a 
horizon. 

2. In the case of SgrA*, the supermassive BH candi¬ 
date at the center of the Galaxy, millimeter and 
infrared observations strongly constrain any possi¬ 
ble thermal blackbody component emitted from the 


surface of this object ( 

Broderick et al. 

2009 |Brod- 

erick & Narayan 

2007 

). Again, the non-detection 


of a thermal component is consistent with the fact 
that the object has no surface but instead possesses 
an event horizon. 

3. In the case of a neutron star, the gas of accretion 
accumulates on the neutron star surface and even¬ 
tually develops a thermonuclear instability. The 
result is a thermonuclear explosion called a Type I 
X-ray burst. The phenomenon seems to be well 
understood and theoretical models agree well with 
observations. These X-ray bursts are observed from 
sources that are supposed to be neutron stars, and 
they have never been observed from BH candi¬ 


dates (Tournear et al. 20031. This should indeed 


be the case if BH candidates have an event hori¬ 
zon, because no gas can accumulate on their sur¬ 


face (Narayan & Heyl 2002 Yuan et al. 2004 1 . 


Bearing in mind the dehnition of a BH, see Subsec¬ 
tion |1.B[ it is fundamentally impossible to test the ex¬ 
istence of an event horizon. These observations may 
instead suggest that BH candidates have an apparent 


horizon. See also Abramowicz et al. (2002) and Bambi 


(2013b) on this point. In particular, there are theoretical 


frameworks in which a BH never forms, instead, only an 


apparent horizon can be created (Bambi et al. 2014a 


2016b Frolov & Vilkovisky 1981). Another issue is that 


the non-detection of a thermal component does not nec¬ 
essarily indicate the absence of a surface. There are sce¬ 


narios, such as those in Dvali & Gomez ( 2013a|b ), in 
which a “classical” BH is actually a Bose-Einstein con¬ 
densate of gravitons. If this were the case, one should 
not expect any blackbody-like component from the BH 
surface. 


II. TESTING THE KERR METRIC 
A. Non-Kerr metrics 

Tests of the Schwarzschild metric in the weak field 
limit are commonly and conveniently discussed within 
the PPN (Parametrized Post-Newtonian) formalism, see 


e.g. Will (2014). We write the most general line element 


for a static and spherically symmetric spacetime in terms 
of an expansion in M/r, namely 


ds^ = - 1- 


2M 




dt^ 


+ (l+ 7 ^ + ...) {dx^ + dy^ + dz^) 


( 9 ) 


where /3 and 7 are two coefficients that parametrize 
our ignorance. In the isotropic coordinates of Eq. ja,. 


the S chwarzschild solution has ,5 = 7=1 (Eddington 
1922[p In Solar System experiments, one assumes that 


/3 and 7 are free parameters to be determined by observa¬ 
tions. Gurrent observational data provide the following 
constraints ([Bertotti et al. 2003[ Williams et al. 2004|) 


1/3-1| < 2.3-10-^ |7-1|<2.3-10 


-5 


( 11 ) 


confirming the validity of the Schwarzschild solution in 
the weak held limit within the precision of current ob¬ 
servations. For the details on how the constraints in 


Eq. (11) are obtained, see the original papers Williams 


et al. (2004) and Bertotti et al. (2003). 


We can attempt to use the same reasoning to test the 
Kerr metric around BH candidates. We hrst consider a 
metric more general than the Kerr solution. Deviations 
from Kerr are parametrized by a number of “deforma¬ 
tion parameters”, which are unknown constants to be 
determined by observations. If the latter require vanish¬ 
ing deformation parameters, the Kerr BH hypothesis is 


The PPN approach is traditionally formulated in isotropic co¬ 
ordinates and the line element of the most general static and 
spherically symmetric distribution of matter is given by Eq. 
When arranged in the more familiar Schwarzschild coordinates, 
the line element becomes 


, 2M , , 2M2 

ds = - (1 —— + iP-'r) 


dV 


-I- + 7?))^ -I- ...^ dr^ + r^dd^ + sin^ ddtf)^ . (10) 
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satisfied. If observations require that at least one defor¬ 
mation parameter is non-vanishing, BH candidates may 
not be the Kerr BHs of general relativity. 

In the Kerr case, it is problematic to find a very gen¬ 
eral solution counterpart of that in Eq. (§. The point is 
that we want to test the strong gravitational field near 
BH candidates and therefore we cannot use an expan¬ 
sion in M/r. Moreover, the metric is now stationary and 
axisymmetric rather than static and spherically symmet¬ 
ric. In the ideal case, we want to have a sufficiently gen¬ 
eral metric that can reduce to any non-Kerr BH in any 
(known and unknown) alternative theory of gravity for 
a specific choice of the values of its deformation param¬ 
eters. At the moment, such a suitable metric with a 
countable number of degrees of freedom does not exist: 
nevertheless there are some proposals in the literature. 

A very popular choice is the Johannsen-Psaltis (JP) 
metric (Johannsen & Psaltis, 2011a). It is an ad hoc 


In (Johannsen & Psaltis 2011a), h is introduced as 


metric and is not a solution of a specific alternative the¬ 
ory of gravity. In Boyer-Lindquist coordinates, the line 
element reads 


ds^ = 


- ( 1 - ^ 1 (1 + /i) dt^ 


AMar sin^ 9 
E 

T,{l + h) 


(1 -|- h) dtdcj) 
dr^ + 


A + ha^ sin^ 9 

, o 2a?Mrs\T? 9 
' + + - 

(E -I- 2Mr) sin^ 9 


sin^ 9d(j)^ , (12) 


where h is 


+ 00 


. / Mr\ /M^Y 

^ ^ f e2fe + e2fc-i-i^w j ( j • (13) 

The metric has an infinite set of deformation parameters 
{ck} and it reduces to the Kerr solution when all the 
deformation parameters vanish. However, eg must vanish 
in order to recover the correct Newtonian limit, while ei 
and £2 are already strongly constrained by experiments 


in the Solar System through Eq. (11) (Cardoso et al. 


2014). The simplest non-trivial metric is thus that with 


63 free and with all the other deformation parameters set 
to zerc 0 


an arbitrary function of r describing deviations from the 
Schwarzschild solution. The rotating JP metric is ob¬ 


tained by applying the Newman-Janis prescription (New¬ 


man et al. 1965 Newman & Janis, 1965), and h becomes 


a function of both r and 9. If we believe that deviations 
from the Kerr geometry must be small, which is consis¬ 
tent with the fact that the known non-Kerr BH solutions 
in alternative theories of gravity are close to the Kerr 
one, h is always a small quantity. In this case, all the 
deformation parameters {cfc} are small and it is possible 
to consider only the unconstrained leading order term, 
neglecting all the others. Alternatively, one may adopt 
a more phenomenological approach, on the basis that 
whatever is not forbidden may be allowed, and consider 
the JP metric with one or a few deformation parameters, 
which are not necessarily small, and simply try to con¬ 
strain that specific choice of the metric. In the following 
sections, I adopt this more relaxed point of view, and 
the deformation parameters will be allowed to have any 
value permitted by observations without the restrictions 
of being small quantities. 

The weak field tests discussed in Subsections lIV.Bl and 
|V.A| can potentially constrain 63 , but can unlikely test 
higher order terms. In particular, if we believe in the in¬ 
terpretation of (Valtonen et al. 2010), 63 would already 


be constrained to a small value. Strong field tests, par¬ 
ticularly in those sources in which the inner edge of the 
disk is very close to the compact object, can instead con¬ 
strain even high order deformation parameters, because 
at the inner edge of the disk M/r can be close to I and 
there is not much difference between lower and higher 
order terms. In the case of strong field tests, with a phe¬ 


nomenological metric like that in Eq. (12), it is common 


to consider just one deformation parameter without the 
restriction of being small. 


® The JP metric is sometimes criticized because of its derivation. 
It is obtained by applying the Newman-Janis prescription, al¬ 
though it is not guaranteed that such an algorithm works outside 
general relativity | |Hansen fc Yun^|2013| . Moreover, the trans¬ 
formation to eliminate some off-diagonal terms is not correct, see 
the discussion in |Azreg-Ainou| ||2011| . 
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An extension of the JP metric is the Cardoso-Pani- 


Rico (CPR) parametrization (Cardoso et al. 2014). The 
line element is 


ds^ = -{l- 


2Mr 


S 

S(l + /in 


A + h'"a'^ sin^ 9 


^ (l + dt^ — 2asin^ 9 v^(l + /i*) (1 + h^) — ^1 — (l + 

2^/{l + h*){l + h’^)- (l- 


dr^ + + sin^ 9 {T. + a^ sin 


dtd(f> 
2Mr 


)(1+ /.*)]} 


#^(14) 


where 


ut ^ f t t Mr\ fM^ 
^ ^ 


+ 00 


E j ’ 


r / r r A/r 

^ = 2_^ \ ^2k + ^2A:+1 

fc=0 




k 


(15) 

(16) 


Now there are two infinite sets of deformation parame¬ 
ters, {el.} and {e^}. The JP metric is recovered when 
el = e} for all k, and the Kerr metric when el = e} = 0 
for all k. 

There are also other proposals in the literature. For 
instance, the quasi-Kerr metric of |Glampedakis fc Babak| 
(2006), which is based on a multipole moment expansion 


and the mass-quadrupole moment of the object is 

Q = Qkbvt: — , 


(17) 


where Qkbtt = is the mass-quadrupole in the 

Kerr case and e is the deformation parameter. Another 


framework is that of the bumpy BHs (Vigeland, 2010 


Vigeland & Hughes, 2010). 


An appealing parametrization, which has not been ex¬ 
plored very far, is represented by the family of regular 


BH metrics of Johannsen (2013b, 2014), which is a gen¬ 


eralization of the family of metrics proposed in [Vigeland] 
et al. ( |2011 ). Here “regular” is to indicate that these 
spacetimes have no naked singularities or closed time-like 
curves outside of the event horizon and the equations of 
motion are separable (as in the Kerr metric, but it is not 
true in general). 

Recently, several others proposed new parametriza- 


tions to test the Kerr metric. Rezzolla & Zhidenko (2014) 


and Konoplya et al. (2016) suggested a parametriza¬ 


tion that seems to fairly well reproduce some known 
non-Kerr BH solutions in alternative theories of gravity 
with a small number of deformation parameters. [Lin etj 
al. (2015) introduced a parametrization suitable for the 


calculations of the electromagnetic spectrum from thin 
accretion disks. Ghasemi-Nodehi & Bambi (2016) dis¬ 


cussed a test-metric in which the free parameters are all 
equal to 1 to recover the Kerr metric and are different 
from 1 in the case of deviations from it. 

Tests of the Kerr metric may also employ the Manko- 


Novikov metric (Manko & Novikov 19921, which was not 


r 


originally proposed to test BH candidates and it is an ex¬ 
act solution of the vacuum Einstein equations with arbi¬ 
trary mass-multipole momentj^ Here the no-hair the¬ 
orem does not apply because the exterior spacetime is 
not regular, namely there are closed time-like curves and 
naked singularities. There are also some variations, such 


as the Manko-Mielke-Sanabria-Gomez solution (Manko 
et al. 2000), which can be extended to include fast- 


rotating objects with J/M^ > 1 (Bambi 2011cI. 

In this subsection we have discussed some model- 
independent parametrizations. An alternative strategy is 
to consider a specific non-Kerr BH solution from a known 
alternative theory of gravity rather than a general ansatz 


like those in (12) and (14). Such an approach is clearly 


much less general. We simply compare two background 
metrics: the Kerr solution against the metric of the al¬ 
ternative theory of gravity of interest. However, there is 
another issue. Usually we know non-rotating solutions in 
alternative theories of gravity, and in a small number of 
cases we know the rotating solution in the slow-rotation 
approximation. Rotating non-Kerr BH solutions in al¬ 
ternative theories of gravity are difficult to find. The 
spin plays a crucial role in the properties of the radiation 
emitted close to BH candidates and therefore, without 
the rotating solution, it is not possible to test the Kerr 
metric, because we are not able to distinguish the effects 
of the spin from those due to possible deviations from 
the Kerr geometry. 


B. Basic properties of non-Kerr metrics 

General statements or properties that hold in the case 
of the Kerr metric may not be true in other backgrounds. 
This may lead to new phenomena with specific astrophys- 
ical implications and observational signatures. 

The uniqueness of the Kerr solution in general relativ¬ 
ity is valid under some specific assumptions and it can 


be violated in the presence of exotic fields (Herdeiro & 


rr 


A particular choice of the mass-multipole of order n fixes the 


value of the current-multipole moment of the same order, which 


is completely determined by the former for given mass-monopole 
and current-dipole moments <|Manko &: Novikov[|1992[|. 












































































11 


Radu 20151. An alternative theory of gravity may po¬ 


tentially have several kinds of BHs, which may be created 
by gravitational collapse from different initial conditions. 
In a similar context, it is possible that astrophysical BH 
candidates are not all of the same type and therefore the 
possible confirmation that a specific object is a Kerr BH 
does not necessarily imply that all BH candidates are 
Kerr BHs. 


1. Horizon 


An event horizon is a null surface in spacetime. If 
we introduce a scalar function / such that at the event 
horizon / = 0, the normal to the event horizon is n'^ = 
d^f and is a null vector. The condition for the surface 
/ = 0 to be null is thus ( Dienerj 2003 Thornburg, 20071 


< 7 ^" ( 5 ^/) ( 9 ./) = 0 . 


(18) 


In general, one can find the event horizon by integrat¬ 


ing null geodesics backward in time, see Diener (2003) 
and Thornburg (2007) for details. In the case of a sta¬ 
tionary and axisymmetric spacetime, the procedure can 
significantly simplify. In a coordinate system adapted to 
the two Killing isometries (stationarity and axisymme- 
try), and such that / is also compatible with the Killing 
isometries, Eq. (fl^ reduces to 


{drfr + 25 "'^ (drf) (def) + {defy = 0 (19) 

in spherical-like coordinates (t, r, 0, (p). The surface must 
be closed and non-singular (namely geodesically com¬ 
plete) in order to be an event horizon and not just a 
null surface. 

If we assume that there is a unique horizon radius for 
any angle 9 (Strahlkorper assumption), we can write / 
as f = r — H{9), where H{d) is a function of 9 only and 
the event horizon is rn = H{9); see Thornburg (2007) 
and [Diener (2003) for details and the limitations of the 


Strahlkorper assumption. The problem is thus reduced 
to finding the solution of the differential equation 


9 


2ff 


re 


dH\ gn 

if) 


dH 

ll9 


= 0 . 


( 20 ) 


The event horizon equation g’"^ = 0 valid in the Kerr 
spacetime in Boyer-Lindquist coordinates only holds 
when the surfaces r = const., which must be closed, 
have a well-defined causal structure, in the sense that the 


surfaces r = const, are null, space-like, or time-like (Jo- 


hannsen 2013a). Such a condition clearly depends on the 


coordinate system, so the event horizon equation g’’’’ = 0 
can only be valid with certain coordinates. 

A Killing horizon is a null hyper-surface on which there 
is a null Killing vector field. In a stationary and ax¬ 
isymmetric spacetime and employing a coordinate sys¬ 
tem adapted to the two Killing isometries, the Killing 


horizon is given by the largest root of 
gttgcjxp ~ gt(f> — 0 ■ 


( 21 ) 


In general relativity, the Hawking rigidity theorem shows 
that the event and the Killing horizons coincide (jHawk- 


mg 


1972), so Eq. (20) and Eq. (21) provide the same 


result. In alternative theories of gravity, this is not guar¬ 


anteed (Johannsen 2013a). 


In general relativity, the event horizon must have 
5^ X R topology, and even this property is regulated by 


mam 


certain theorems (Hawking, 1972 Jacobson & Venkatara- 


1995). For instance, toroidal horizons can form. 


but they can only exist for a short time, in agreement 


with these theorems (Hughes et ah, 1994). If we want 


to test the Kerr metric and general relativity, we can¬ 
not exclude scenarios with BHs with a topologically 


hannsen 2013a). 


non-trivial event horizon (Bambi & Modesto 2011 Jo- 


2. Particle orbits 


In a generic stationary and axisymmetric spacetime, 
there are three constants of motion, namely the mass m, 
the energy E, and the axial component of the angular mo¬ 
mentum Lz. The Kerr solution is a Petrov type D space- 
time and therefo re there is a fo urth constant Q called the 
Carter constant (Carter 1968p If the spacetime has a 
forth constant of motion, it is always possible to choose 
a coordinate system in which the equations of motion 
are separable. In the Kerr metric, this is the case for 
the Boyer-Lindquist coordinates and the calculations of 
the photon trajectories from the observer to the region 
around the BH can be reduced to the calculations of some 
elliptic integrals, which is computationally more efficient. 
In generic spacetimes with a Carter-like constant, these 
calculations are already more complicated and, in the 
simplest case, the elliptic integrals of the Kerr metric be¬ 


come hyper-elliptic integrals (Enolski et al. 2011). If the 


spacetime has no Carter-like constant, it is necessary to 
solve a system of coupled second order differential equa¬ 
tions. 

In accretion disk models and astrophysical measure¬ 
ments, the ISCO radius plays an important role. For 


instance, in the Novikov-Thorne model (Novikov & 


Thorne 1973), which is the standard set-up for thin ac¬ 


cretion disks, the inner edge of the disk is at the ISCO 
radius and it is the crucial parameter to measuring the 
spin or possible deviations from the Kerr solution. In the 
Kerr metric, equatorial circular orbits are always ver¬ 
tically stable, whereas they are only radially stable for 


In the Schwarzschild limit (a* = 0), the Carter constant reduces 
to Q = — Lz, where L is the total angular momentum. In the 

Kerr metric, there is not a direct physical interpretation of Q. 
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T > nsco- In non-Kerr metrics, the picture is more 
complicated. The ISCO may be either radially or ver- 

| 2011 a|b| ), 


tically marginally stable (Bambi & Barausse 


and it is also possible to have disconnected stable regions 


among which equatorial orbits are unstable (Bambi & 


Lukes-Gerakopoulos 2013 1 . Moreover, in the Kerr met¬ 
ric the usual picture is that a particlj^ of the accretion 
disk reaches the ISCO and then quickly plunges onto the 
BH, without emitting much radiation after leaving the 
ISCO. In non-Kerr metrics, it is possible that the parti¬ 
cle plunges to a region inside the ISCO, but it cannot be 


immediately swallowed by the BH (Bambi & Barausse 


2011b I. On the contrary, there may be an accumulation 
of gas between the ISCO and the compact object, and 
the gas has to lose additional energy and angular mo¬ 
mentum before plunging to the central body. There are 
also spacetimes in which there is no ISCO, namely the 
orbits are always stable. For instance, this is the case of 
some metrics describing the spacetime around certain ex¬ 


otic non-interacting matter (Bambi & Malafarina 2013 


Kesden et al. 

2005 

Macedo et al. 

2013 


In the Kerr metric, the ISCO radius is located at the 
minimum of the energy of equatorial circular orbit At 
larger radii, the specific energy monotonically increases 
to approach 1 at infinity. From the ISCO radius to 
smaller radii, the specific energy monotonically increases 
to diverge to infinity at the photon orbit. At the radius 
of the marginally bound circular orbit, which is located 
between the photon radius and the ISCO radius, the spe¬ 
cific energy of a test-particle is if = 1. When the ISCO 
is marginally vertically unstable, the energy of equato¬ 
rial circular orbits may be a monotonic function without 
minimum (the energy decreases as the radial coordinate 
decreases). In this case, there is no marginally bound 
circular orbit. Accretion disk structures around similar 


BHs are qualitatively different (Li & Bambi 2013a). 


Since astrophysical observations are often sensitive to 
the position of the ISCO, it is useful to have an idea 
of the correlation between the spin and the deformation 


damental frequencies in the JP metric are correlated to 
the position of the ISCO radius. The left panel of Fig. 
shows the contour levels for risco = 2 M, 3 M, 4 M, 5 M, 
and 6 M on the plane spin parameter a* and JP defor¬ 
mation parameter As already mentioned, throughout 
this article I do not impose that the deformation param¬ 
eter must be a small quantity. 

For very large and positive £ 3 , the ISCO is at a larger 
radius because of the strong instability along the verti¬ 
cal direction. The contour levels of the ISCO can give a 
simple idea of which spacetimes may look similar in as¬ 
trophysical observations. The right panel in Fig. shows 
the Novikov-Thorne radiative efficiency ryNx = 1 ~A'iscOj 
where Aisco is the specific energy of a test-particle at 
the ISCO radius. The contour levels are for t^nt = 0.05, 
0.10, 0.15, 0.20, 0.25, and 0.30. ? 7 nt is indeed a more 
appropriate quantity than the ISCO radius. First, ryNx 
is independent of the coordinate system (the radial coor¬ 
dinate of the ISCO has no physical meaning) and it can 
be directly measured. Second, it is indeed a better esti¬ 
mator to figure out which spacetimes may look similar in 
astrophysical observations (see the two panels in Fig. 
and compare the shape of their contour levels with, e.g., 
the shape of the constraints in Figs. and [^. 

Finally, note that non-Kerr metrics often have some 
pathological features for some choices of the deforma¬ 
tion parameters. Naked singularities, regions with closed 
time-like curves, etc, are possible. Some caution has to 
be taken. However, even pathological spacetimes can be 
used to test the Kerr metric if we assume that the space- 
time solution is only valid outside of some interior region. 
For example, the interior would be different because of 
matter source terms. In the case of a compact object 
made of exotic matter, the vacuum solution would hold 
up to the surface of the object, while at smaller radii 
the metric would be described by an interior solution. 
Pathological features like naked singularities may also be 


removed by unknown quantum gravity effects (Cimon & 


Hofava, 2009). A study of the pathologies in these space¬ 


Johannsen & Psaltis 

(20131 showed how iron times is reported m 

Johannsen ( 

2013a 


lines in the JP metric with different o» and £3 but the 
same ISCO radius are very similar, as well as orbital fun- 


It is worth pointing out that, in this discussion and in what 
follows, “particle” is used to refer to a parcel of gas. 

This can be seen as follows. The specific energy of a test-particle 
in an equatorial circular orbit is given, for instance, in Eq. (2.12) 
in |Bardeen et alT] | |1972[ |. Its derivative with respect to the radial 
coordinate is 


dE 


- GMr + - ‘ia? 


dr 2r7/4 (^3/2 _ 


( 22 ) 


The minimum can be found from dEjdr = 0, which is the same 
equation as that for the orbital stability, see e.g. [Bardeen et al.| 
l|1972| or |Chandrasekhar| ||1985| . 


C. Violation of the Kerr bound |a*| < 1 

A fundamental limit for a Kerr BH is the bound 
|a*| < 1. As seen from Eq. 0 , for I a* I > I there is 
no horizon, and the Kerr metric describes the spacetime 
of a naked singularity. According to the cosmic censor¬ 
ship conjecture, naked singularities cannot be created by 


gravitational collapse (Penrose 1969), even if we know 


some counterexamples in which a naked singularity ex¬ 


ists for an infinitesimal time (Joshi & Malafarina 2011 1 . 


The question is thus if BH candidates may be objects 
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FIG. 3 Contour levels of the ISCO radius risco (left panel) and of the Novikov-Thorne radiative efficiency ?7 nt = 1 — l^isco 
(right panel) in the JP metric with non-vanishing deformation parameter £3. The black dotted line separates spacetimes with 
a regular event horizon (on the left of the line) from those with topologically non-trivial event horizons and naked singularities 
(on the right). See the text for more details. 


with |a*| > 

First, an object with |a*| ^ 1 is a fast-rotating body 
only in the case where it is very compact. For instance, 
the spin parameter of Earth is about 10^ and there is no 
violation of any principle because the vacuum solution 
holds up to that of the Earth surface. The Earth’s ra¬ 
dius is rEarth ~ 6,400 km, which is much larger than its 
gravitational radius = M « 4.4 mm. 

Second, the Kerr metric with |a*| > 1 is not a viable 
astrophysical scenario. For instance, [Giacomazzo et al.| 
(2011) found that in general relativity stellar models with 


I a* I > 1 do not collapse without losing angular momen¬ 
tum. If we consider an existing Kerr BH and we try to 
overspin it up to |a*| > 1, we fail (Barausse et al., 2010). 


The same negative result is found by considering the col¬ 
lision of two BHs at the speed of light in full non-linear 
general relativity (Sperhake et al. 2009). None of these 


studies proves that it is impossible to create a Kerr space- 
time with I a* I > 1, but at least they suggest that naked 
Kerr singularities may not be physical merely because 
they may not be created in an astrophysical process. As¬ 
suming it is somehow possible to create a similar object. 


Pani et al. (2010) showed that the spacetime would be 


unstable independently of the boundary conditions im¬ 
posed at the excision radius r, which means that possible 
unknown quantum gravity corrections at the singularity 
cannot change these conclusions. Very compact objects 
described by the Kerr solution with |a*| > 1 are thus 
unlikely astrophysical viable BH candidates. 

Third, if the measurement of the spin parameter of 
a BH candidate gave a value larger than 1 assuming 


the Kerr metric, the measurement would be presumably 
wrong, but it would be a clear indication of new physics. 
As previously discussed, the Kerr metric with |a*| > 1 is 
not a viable astrophysical scenario. The metric around 
the BH candidate should thus be different from that of 
the Kerr metric. Since spin measurements strongly de¬ 
pend on the exact background metric, the measurement 
would be wrong and it is possible that the actual value 
of the spin parameter is instead smaller than 1. 

Last, note that |a*| = 1 is a critical value only in the 
Kerr metric. In the case of non-Kerr BHs, the critical 
bound is typically different, and it may be either larger 
or smaller than 1, depending on the spacetime geometry. 
Moreover, there are examples in which one can create a 


non-Kerr compact object with |a*| > 1 (Bambi 2011c| d) 
or in which one can overspin a BH up and “destroy”|^° 


its horizon (Li & Bambi, 2013b). 


D. Important remarks 

In order to test the Kerr metric with electromagnetic 
radiation, we need to study the properties of the radia¬ 
tion emitted by the gas in the accretion disk or by stars 
orbiting the BH candidate. Assuming we are in a metric 


theory of gravity (Will 2014), namely that test-particles 


follow the geodesics of the spacetime, the metric provides 
all the answers. The spectrum of the BH candidate de¬ 
pends on the motion of the gas in the accretion flow or 


We remind the reader that here, as well as in the remainder of 
the article, a* = J/M^ is the dimensionless spin parameter of 
the compact object, not merely a parameter of the metric. 


An event horizon and a BH cannot be destroyed by definition. 
Here we start from a stationary spacetime describing a BH and 
we introduce some small particles to destroy the BH of the sta¬ 
tionary metric. The small particles thus completely change the 
causal structure of the spacetime. 
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of orbiting stars and by the propagation of the photons 
from the point of emission to the distant observer. We 
thus note the following: 

1. With this approach, we test the Kerr metric, just 
like with the PPN formalism we can test the 
Schwarzschild solution. We do not directly test the 
Einstein equations. For instance, we cannot distin¬ 
guish a Kerr BH of general relativity from a Kerr 
BH in an alternative metric theory of gravity, be¬ 
cause the motion of particles and photons is the 
same. Actually the Kerr metric is quite a com¬ 


Within the corona-disk model with lamppost geometry, 
the set-up is shown in Fig. ([Martocchia fc Matt 1996 


mon solution to many gravity theories (Psaltis et 


al. 20081. However, the observation of a non-Kerr 


BH could rule out the Einstein equations because, 
within general relativity, astrophysical BHs should 
be well described by the Kerr solution. 

2. If we want to directly test the Einstein equa¬ 
tions, we should use gravitational (rather than 
electromagnetic) waves, because their emission is 
governed by the Einstein equations ( ]Barausse fc 
Sotiri^ 2008| ) . 


3. Even if we assume the Kerr metric, analyze the 
data, and obtain a very good fit, it is not enough 
to claim that the object is a Kerr BH. This is made 
clearer in the next sections. There is typically a de¬ 
generacy among the parameters of the model, and 
in particular between the spin and possible devia¬ 
tions from the Kerr background. For this reason 
we adopt a PPN-like approach. 


III. X-RAY OBSERVATIONS 

The X-ray radiation in the spectrum of BH candidates 
is thought to be generated in the vicinity of these ob¬ 
jects. Some features depend on the strong gravitational 
field near the BH candidate and therefore, if properly 
understood, they may be used to test the Kerr BH hy¬ 
pothesis. 

There is some difference between stellar-mass and 
supermassive BH candidates, because of the different 
masses and environments. Generally speaking, BH can¬ 
didates are observed in different spectral states, which 
correspond to different geometries and emission mecha¬ 


nisms of their accretion flow (Belloni, 20101. In the case 


of stellar-mass BH candidates, a source can change its 
spectral state on a timescale of weeks or months. For 
supermassive BH candidates, the timescales are too long 
and the source can only be observed in its current spec¬ 
tral state. Note that the spectral state classification is 
still a work in progress, some spectral states and their 
physical interpretation are not yet well understood, and 
different authors may use a different nomenclature. The 
key-point in any measurement is to have the correct as¬ 
trophysical model. 


Matt et al. 1991). The accretion disk is geometrically 
thin and optically thick. It radiates like a blackbody 
locally and as a multi-color blackbody when integrated 
radially. The temperature of the disk depends on the BH 


mass and the mass accretion rate [see e.g. Bambi et al. 
(2016a)], and the disk’s thermal spectrum is in the X-ray 
band only for stellar-mass BH candidates. The “corona” 
is a hotter, usually optically-thin, electron cloud which 
enshrouds the central disk and acts as a source of X-rays, 
due to inverse Compton scattering of the thermal photons 
from the accretion disk off the electrons in the corona. 
The corona is often approximated as a point source lo¬ 
cated on the axis of the accretion disk and just above the 
BH. This arrangement is often referred to as a lamppost 
geometry. The lamppost geometry requires a plasma of 
electrons very close to the BH and such a set-up may be 
realized in the case of the base of a jet. However, other 
geometries are possible, and an example is the family of 
“sandwich models” (Schnittman & Krolik, 2010). The 


direct radiation from the hot corona produces a power- 
law component in the X-ray spectrum. The corona illu¬ 
minates also the disk, producing a reflected component 
and some emission lines, the most prominent of which is 
usually the iron Ka line. 

X-ray techniques to test the Kerr metric are dis¬ 
cussed in the next subsections: continuum-fitting method 


(HI.A), iron Ka line with time-integrated (HI.B) and 
time-resolved data (HI.C), quasi-periodic oscillations 
(QPOs, [lII.D ), and polarization measurements (HI.E). 

Currently, the two leading techniques to probe the 
spacetime geometry around BH candidates with X-ray 
measurements are the study of the thermal spectrum of 


thin disks (continuum-fitting method) (Li et 

al. 

2005 

McClintock et al. |2014| Zhang et al.| |1997 

) and the 

analysis of the iron Ka line (Brenneman & Reynolds 

2006 

Fabian et al. 1989 

Reynolds 

2014 

). Both tech- 


niques have been developed to measure the spin param¬ 
eter of BH candidates under the assumption of the Kerr 
background, but they can be naturally extended to test 
the Kerr metric. Current spin measurements of stellar- 
mass BH candidates reported in the literature under the 
assumption of the Kerr background are summarized in 
Tab. |T| Iron line spin measurements of supermassive BH 
candidates under the assumption of the Kerr background 
are reported in Tab. 

QPOs are seen as peaks in the X-ray power density 
spectra of BH candidates and they may be used to mea¬ 
sure the properties of the metric around these objects. 
For the time being, we do not know the exact mechanism 
responsible for these oscillations, and therefore QPOs 
cannot yet be used to test fundamental physics. Different 
models provide different results. However, QPOs are a 
promising tool for the future, because the value of their 
central frequency can be measured with high precision. 
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Power-Law Component 



FIG. 4 Corona-disk model with lamppost geometry. The BH 
is surrounded by a geometrically thin and optically thick ac¬ 
cretion disk. The accretion disk radiates like a blackbody lo¬ 
cally and as a multi-color blackbody when integrated radially. 
The hot electron cloud called corona acts as an X-ray source 
and it is located just above the BH. The power-law compo¬ 
nent represents the direct radiation from the hot corona. The 
latter illuminates also the disk, producing a reflection compo¬ 
nent and some emission lines, the most prominent of which is 
usually the iron Ka line. Figures courtesy of Jiachen Jiang. 


The observation of the X-ray polarization of the ther¬ 
mal spectrum of thin accretion disks is another poten¬ 
tial technique to test the Kerr metric. There are cur¬ 
rently no space missions with polarization detectors, but 
there are a few proposals that may be operative within 


10 years: the Chinese-European mission eXTP (Zhang 
et al. 2016[), the European project XIPE (Sofhtta et al. 


2013 ), or the two NASA missions IXPE (jWeisskopf et al. 
2008| and PRAXY^ 


A. Continuum-fitting method 


The continuum-fitting method consists of the analysis 
of the thermal spectrum of geometrically thin and opti- 


cally thick accretion disks 1 

Li et al. 

2005 

McClintock 

et al. 

to 

o 

Zhang et al. 

1997 

|. The standard theoreti- 


Thorne 1973 Page & Thorne 19741. The plane of the 
disk is assumed to be perpendicular to the BH spiif^ 
The particles of the gas move on nearly geodesic circular 
orbits, and the inner edge of the disk is at the ISCO ra¬ 
dius. The latter assumption plays a crucial role and it is 


http://asd.gsfc.nasa.gov/praxys/ 

If the BH formed from the supernova explosion of a heavy star 
in a binary, its spin would more likely already be aligned with 
the orbital angular momentum vector of the system (at least 
in the case of a symmetric explosion without strong shock and 
kick) | |Fragos et al.[[2010[ |. If this is not the case, the inner part 
of the accretion disk can still be expected to be on the plane 
perpendicular to the BH spin as a consequence of the Bardeen- 
Pette rson effect jBardeen &: P etterson| |1975| [Kumar &: Pringle| 

|1985| [Steiner & McClintockj |2012|l. 


confirmed by some observations that show that the inner 
edge of the disk does not change appreciably over several 
years when the source is in the thermal state. The most 
compelling evidence comes from LMC X-3. The analy¬ 
sis of many spectra collected during eight X-ray missions 
and spanning 26 years shows that the radius of the inner 
edge of the disk is quite constant (Steiner et al. 20101. 
The most natural interpretation is that the inner edge is 
associated with some intrinsic property of the geometry 
of the spacetime, namely the radius of the ISCO, and it 
is not affected by variable phenomena like the accretion 
process. 

An accretion disk meets these conditions when the 
source is in the so-called high/soft state. Here “high” 
refers to the accretion luminosity, which must be higher 
than 5% of the Eddington limit. “Soft” refers to the 
fact that the soft X-ray component, corresponding to the 
thermal spectrum of the disk, is the dominant one. This 
is true only for stellar-mass BH candidates, because the 
temperature of the disk depends on the BH mass and 
the mass accretion rate. As it turns out, the thermal 
spectrum of a thin disk is in the X-ray range for stellar- 
mass BH candidates and in the UV/optical bands for 
supermassive BH candidates. In the latter case, extinc¬ 
tion and dust absorption limit the ability to make ac¬ 
curate measurements. The continuum-fitting method is 
thus normally applied to stellar-mass BH candidatej^ 
Thin disks are present when the accretion luminosity is 
5-30% of the Eddington limit (McClintock et al. 2006|). 
For the validity of the method, see e.g. McClintock et al!] 
(2014) and references therein. 

In the Kerr background, the thermal spectrum of a 
thin disk only depends on 5 parameters: the BH mass 
M, the mass accretion rate M, the inclination angle of 
the disk with respect to the line of sight of the observer 
i, the distance d of the source, and the spin parameter 
of the BH a*. The impact of the model parameters on 
the shape of the spectrum is shown in Fig. If it is 
possible to have independent measurements of M, i, and 
d, one can ht the thermal component of the spectrum and 
determine the BH spin a* and the mass accretion rate M. 
Current spin measurements with the continuum-fitting 
method are reported in Tab. |l] Some of these objects are 
not dynamically confirmed BH candidates, so their mass 
is not estimated from the motion of the companion star. 

The technique can be naturally extended to non-Kerr 
backgrounds. Preliminary studies considered some spe¬ 
cific non-Kerr objects and did not take into account all 


the relativistic effects (Bambi & Barausse 

2011 a 

Harko 

et al. 

2009a|b 

2010 

Pun et al. 

2008a|b 

Torres 

2002 ). 


The continuum-fitting method was applied to supermassive BH 
candidates in some very special cases in [Czerny et al.[ | |2011[ ) and 
[Done et al.[l|2013^. 













































































































16 


BH Binary o* (CF) a* (iron) Principal References 


GRS 1915-105 > 0.98 0.98 ±0.01 

McClintock et al. (2 

OOol; 

Miller et al. 

|2013j) 

Cyg X-1 > 0.98 0.971° °^^ 

Fabian et al 

( 2012 ) 

; Gou 

et al. 

(2011 

2014) 

LMC X-1 0.92 ± 0.06 0.97lo “5 


Gou et al. 

(20091 

Steiner et a 

il ( 2012 b| 

GX 339-4 < 0.9 0.951^1 Gardaetal. (2 

015 i; K 

olehn 

lainen 

& Done 

(20101 

MAXI J1836-194 — 0.88 ± 0.03 

M33 X-7 0.84 ± 0.05 — 

4U 1543-47 0.80 ±0.10* — 

Swift J1753.5 — 0.76l2;(g 

XTE J1650-500 — 0.84 ~ 0.98 

IC 10 X-1 > 0.7 — 


Reis e 

t al. ( 

2012 ) 



L 

in et al 

(200 

8 201 c 

1 


Shafee et al. 

(2006) 


Reis et al. ( 

2009) 



Walton et al. 

(2012 


Steiner et al. 

2016 

GRO J1655-40 0.70 ±0.10* >0.9 

Reis et a 

.(2009 

1 ; Sha 

Fee et 

il. (2006 

1 

GS 1124-683 0.63t[|;(^ — 

XTEJ1752-223 — 0.52 ±0.11 

XTE J1652-453 — < 0.5 

XTE J1550-564 0.34 ± 0.28 0.55lg;22 

LMC X-3 0.25 ±0.15 — 

H1743-322 0.2 ± 0.3 — 

A0620-00 0.12 ±0.19 — 

XMMU J004243.6 < -0.2 — 


(Chen et al. 

20161 

1 

Reis et al. ( 

2011 ) 


Chiang et al. 

(2012 


Steiner et al. 

(2011 


Steiner et al. 

(2014 

c 

jteiner et al. 

2012 a 


Gou et al. ( 

20101 



M 

iddleton et a 

1 . (201 

4) 


TABLE I Summary of the continuum-fitting (CF) and iron line measurements of the spin parameter of stellar-mass BH 
candidates under the assumption of the Kerr background. See the references in the last column for more details. *These 
sources were studied in an early work on the continuum-fitting method, within a more simple model, so the published 1-a error 
estimates are doubled following McClintock et al. (20141. 


A ray-tracing code that includes all the relativistic ef¬ 
fects and can compute the thermal spectrum of a thin 
disk in a generic stationary and axisymmetric spacetime 
was presented in Bambi] (2012a|. 

The impact of a possible non-vanishing deformation 
parameter on the thermal spectrum of a thin disk is 
shown in the bottom right panel in Fig. Even with¬ 
out a quantitative analysis, it is clear that the effect of 
the spin and of the deformation parameter is very sim¬ 
ilar. The point is that the cut-off of the spectrum is 
determined by the position of the inner edge of the disk, 
which is assumed to be at the ISCO radius, or, more pre¬ 
cisely, by the radiative efficiency in the Novikov-Thorne 
model (Kong et al. 2014). The radiative efficiency in the 
Novikov-Thorne model is 


ryNT = 1 — Aisco , (23) 

where Eisco is the specific energy of a test-particle at the 
ISCO radius. In the Kerr metric, 77 nt ranges from 0.057 
for a non-rotating BH (a* = 0) to about 0.42 for a maxi¬ 
mally rotating BH and a corotating disk (a* = 1). In the 
case of retrograde disks, decreases as the spin param¬ 
eter increases until about 0.038 for a* = 1. In the Kerr 
metric, there is a one-to-one correspondence between 77 NT 
and a*, and therefore the measurement of the former can 
provide an estimate of the BH spin parameter. If we relax 
the Kerr BH hypothesis and we allow for a non-vanishing 
deformation parameter, the same value of the radiative 


efficiency can be obtained for different combinations of 
the spin parameter and the deformation parameter. The 
result is that there is a degeneracy and it is impossible to 
measure both the spin and the deformation parameter. 
In general, we can measure only a combination of them. 

In Kong et al. (2014), we reconsidered the spin mea¬ 
surements with the continuum-fitting method reported 
in the literature and under the assumption of the Kerr 
background. We obtained the constraints on the spin 
parameter - deformation parameter plane within the JP 
background. Some examples are shown in Figs. and 
0 In the case of A0620-00, LMC X-3, XTE JI550-564, 
and 4U 1543-47, there is a specihc measurement in the 
Kerr metric, and this can be translated into an allowed 
region on the spin parameter - deformation parameter 
plane. These constraints are obtained with the following 
method that can be justihed a posteriori [more details 
can be found in the original paper Kong et al. (2014)]. 
Instead of working on some observational data, we adopt 
the theoretical spectrum of a disk around a Kerr black 
hole with the parameter values reported in the litera¬ 
ture (M, i, and d inferred without the assumption of the 
Kerr metric, and a* and M obtained with the continuum¬ 
fitting method and assuming the Kerr metric). Such a 
spectrum is then compared with the theoretical spectra 
computed in spacetimes in which a*, M, and €3 vary. 
Employing a analysis, the red solid lines in Fig. are 
the best fits for a fixed £ 3 . The blue dashed lines are the 
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AGN _ a, (iron) I/Boi/7/Edd 


IRAS 13224-3809 
Mrk no 
NGC 4051 
MGG-6-30-15 
IH 0707-495 

NGC 3783 
RBS 1124 
NGC 1365 
Swift J0501.9-3239 
Ark 564 
3C 120 
Ark 120 

Ton SI80 
IH 0419-577 
Mrk 509 

IRAS 00521-7054 
3C 382 
Mrk 335 
Mrk 79 
Mrk 359 
NGC 7469 
Swift J2127.4+5654 
Mrk 1018 
Mrk 841 
Fairall 9 


> 0.995 

0.71 

> 0.99 

0.16 ±0.04 

> 0.99 

0.03 

> 0.98 

0.40 ±0.13 

> 0.98 

~ 1 

> 0.98 

0.06 ± 0.01 

> 0.98 

0.15 

n 07 + 0.01 
-0.04 

0 06+° °® 
U.UD_o 04 

> 0.96 

— 

'J-yO-0.06 

> 0.11 

> 0.95 

0.31 ±0.20 

0.94 ±0.01 

0.04 ±0.01 

n 01 + 0.02 
'^•^J^-0.09 

Q 1 +3.2 
^•-*-- 1.6 

> 0.88 

1.3 ± 0.4 

0 . 86 t°:“ 

— 

> 0.84 

— 

'^■•'^-0.04, 

— 

0 70+012 
u- ‘U-O.Ol 

0.25 ±0.07 

0.7 ±0.1 

0.05 ±0.01 

0 7+00 
o’-'-o.s 

0.25 

0.69 ±0.09 

— 

0.6 ± 0.2 

0.18 ±0.03 

0 O'*'® '* 

U.O -0 7 

0.01 

> 0.56 

0.44 

0.52tVil 

0.05 ±0.01 


Principal References 


Brenneman et al. 


(2011I; 



Walton et al. 

(2013 

1 

Walton et al. 

(2013 




fiatrick et al. ( 

2011 a 

1 

Brenneman & R 

eynolds 

(2006 

; Min 

iutti et al. 

(2007 

1 

de La Galle Per 

3 z et al. 

(2010 

j Wal 

ton et al. 

2013 

7 


Zoghbi et al. 

2010 





Patrick et al. 


Walton et al. 


Brenneman et al. 


120131; 


(2013 I 


( 2011 al 


Risaliti et al. 


Walton et al. 


Walton et al. 


Lohfink et al. 


Nardini et al. 


Walton et al. 


Walton et al. 


( 20111 ; 


(2013 I 


(2013 


(20131 


(2013 


20091 


Patrick et al. 


Walton et al. 


Walton et al. 


Tan et al. 


(20131 


(20131 


(20131 


(20131 


(20121 


Walton et al. 


( 2011 bl; 


(20131 


Patrick e 


il. 

( 2011 b 1 

W 

^al 1 

on 

et 

al. 

(2013 


G 

alio et al. 

20 

05 

20 

11 

1 


Walton et al. 


Patrick et al. 


Miniutti et al. 


(20131 


( 2011 b| 


We 

ilton 

J. 

it 

al. 

(2013 

We 

ilton 

3t 

al. 

(2013 


Lohfink et al. 


( 20121 ; 


Schmoll et al. 


(2011b I 


Patrick et al. 


j 2011 b(; 


(20091; Walton et al. (20131 


TABLE II Summary of the iron line measurements of the spin parameter of supermassive BH candidates under the assumption 
of the Kerr background. See the references in the last column and [Brenneman |( (M^ for more details. 


boundaries of the allowed regions at l-tr. However, in 
this way we are assuming that these sources have vanish¬ 
ing 63 . For a fixed es, the l-cr error on the spin is given 
by the dash-dotted green lines, which should better ap¬ 
proximate the actual l-cr error if we analyzed the real 
data. This approach can be justified a posteriori because 
the constraints provided by the dash-dotted green lines 
and by the dashed blue lines are very similar consider¬ 
ing the spin uncertainty. There is a quasi-degeneracy in 
the theoretical prediction of the spectra and therefore a 
spin measurement inferred in the Kerr background can be 
translated into an allowed region on the spin parameter 
- deformation parameter plane. 

Fig.j^shows the constraints on €3 from GRS 1915-1-105 
and Cygnus X-1. In both cases, the Kerr measurement 
is a* > 0.98, and therefore in the JP metric we have 
an allowed region in which the spectra look more like a 
Kerr BH with a* > 0.98 and an excluded region in which 
the spectra are more like that of a Kerr BH with a* < 
0.98. The difference between the two allowed regions is 


only due to the different inclination angle, i = 66° for 
GRS 1915-1-105 and i = 27.1° for Gygnus X-1. Note the 
similarity of the shapes of the allowed regions in Figs. 
and with the contour levels of t^nt in the right panel in 

Fig-i 


If we have a BH candidate that looks like a very fast- 
rotating Kerr BH, similar to that in Gygnus X-1 in the 
example in Fig. it is sometimes possible to constrain 
the deformation parameter. Another example with a dif¬ 
ferent background metric was shown in Bambi (2014a I. 
The reason is that, in general, if one considers very large 
deviations from Kerr, in both directions in the deforma¬ 
tion parameter, the ISGO radius increases and therefore 
the Novikov-Thorne radiative efficiency decreases. The 
result is that very deformed objects cannot mimic a fast- 
rotating Kerr BH. However, this is not a general state¬ 
ment, and some deformations may be extremely large. 
An example of the latter case is the CPR deformation pa¬ 
rameter Cg, which cannot be constrained with Gygnus X- 
1 ( ]Bambi 2014bI, see Fig. 
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^obs 


^obs 



^®9 ^obs 



^^9 ^obs 



^®9 ^obs 



^®9 ^obs 


FIG. 5 Impact of the model parameters on the thermal spectrum of a thin disk: mass M (top left panel), mass accretion rate 
M (top right panel), viewing angle i (central left panel), distance d (central right panel), spin parameter a* (bottom left panel), 
and JP deformation parameter £3 (bottom right panel). When not shown, the values of the parameters are: M — 10 Mq, 
M = 2 ■ g s“^, d = 10 kpc, i = 45°, a* = 0.7, and £3 = 0. M in units of Mq, M in units of 10^® g s“^, d in kpc, flux 
density in photons keV“^ cm“^ s“^, and photon energy Eohs in keV. 


In principle, one could find a source with a thermal 
spectrum harder than that which is expected for a Kerr 
BH with a* = 1. This would essentially correspond to 
a spacetime in which the Novikov-Thorne radiative effi¬ 
ciency exceeds the Kerr BH bound 77 NT = 0.42 and could 


be an indication of deviations from the Kerr geometry. 
For the time being, there are no observations of this kind 
and therefore all the data are consistent with the Kerr 
metric. 

The continuum-fitting method is probably the most 
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A0620-00 


LMC X-3 




XTE J1550-564 4U 1543-47 




FIG. 6 Continuum-fitting constraints on the JP parameter es from A0620-00 (top left panel), LMC X-3 (top right panel), 
XTE J1550-564 (bottom left panel), and 4U 1543-47 (bottom right panel). The spacetimes along the red solid lines cannot 
be distinguished and represent the central values of the measurement. In every panel, the dash-dotted green lines are the 
boundaries of the allowed region (l-cr error) along the red line inferred within the analysis of |Kong et al.| ( [2014[ ) from current 
X-ray measurements. The blue dashed lines are the same boundaries inferred if the best lit were exactly for £3 = 0. The dotted 
black curve on the right of each panel separates spa cetimes with a regu lar e xterior (on the left of the curve) from those with 
naked singularities (on the right of the curve). From Kong et al. ( 2014|. See Kong et al. (20141 and the text for more details. 


robust technique among those available today. Its as¬ 
sumptions have been tested in a number of studies; for 
instance, the observed temporal constancy of the accre¬ 
tion disk’s inner radius in the thermal state supports the 


assumption that the inner edge is at the ISCO (Steiner 


et al., 2010). However, the method also has some weak 


points. It cannot be applied to AGN, which represent the 


majority of the BH candidates, because their disk tem¬ 
perature is in the UV/optical bands. The measurements 
of M, i, and d from optical observations are sometimes 
difficult and may be affected by systematic effects. Cor¬ 
rections for non-blackbody effects are usually taken into 
account by introducing the color factor, which is obtained 
from disk atmosphere models and there is not a unani- 
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FIG. 7 Continuum-fitting constraints on the JP parameter es from GRS 1915+105 (left panel) and Cygnus X-1 (right panel). 
The allowed regions are those inside the dashed blue lines. The dotted black curve separates spacetime with a regular exterior 
(on the left of the curve) from those with naked singularities (on the right of the curve). From liTong et ah ( 2014| ). See Kong 
|et al.| (|2014[) and the text for more details. 


mous consensus on their reliability flDavis et al. 2005 


Davis & Hubeny 20061. 


Assuming the systematics are under control, the ther¬ 
mal spectrum of a thin disk has a very simple shape, 
and it cannot provide much information on the space- 
time geometry around the BH candidate. If we assume 
the Kerr metric, we can determine the spin parameter. If 
we have just one possible non-vanishing deformation pa¬ 
rameter, we meet a degeneracy and, in general, we can¬ 
not constrain the spin and possible deviations from the 
Kerr solution at the same time. If we have a source that 
looks like a very fast-rotating Kerr BH, we can constrain 
some deformation parameters (e.g. the JP deformation 
parameters ek), but other deviations from the Kerr ge¬ 
ometry cannot be constrained (e.g. the CPR deforma¬ 
tion parameters ep. The reason is that the spectrum is 
simply a multi-color blackbody spectrum without addi¬ 
tional features. Different parameters of the model have 
a quite similar impact on the shape of the spectrum and 
therefore there is a strong parameter degeneracy. The 
best that we can do is to combine the continuum-fitting 
measurements with other observations to break the pa¬ 
rameter degeneracy. 


Cygnus X-1 



a. 


FIG. 8 Gontinuum-fitting constraints on the CPR parameter 
£3 from Cygnus X-1. The allowed region is between the two 
dashed blue lines. See the text and the original paper |Bambi| 
(2014b I for more details. 


B. Iron line spectroscopy 

The illumination of a cold disk by a hot corona pro¬ 
duces a reflection component as well as some spectral 
lines by fluorescence in the X-ray spectrum of the source. 


The most prominent line is usually the iron Ka line, 
which is at 6.4 keV in the case of neutral atoms and 
it shifts up to 6.97 keV for H-like iron ions. This line is 
intrinsically narrow in frequency, while the one observed 
in the X-ray spectrum of BH candidates appears broad¬ 
ened and skewed. The interpretation is that the line is 
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FIG. 9 Photon count (in black) and fit without the iron 


line (in red) of MGC-6-30-15. From Brenneman & Reynolds 


(20061. ©AAS. Reproduced with permission. 


strongly altered by relativistic effects, which produce a 
very characteristic profile first studied in [Fabian et al.| 
(1989). This interpretation is currently well supported 
by reverberation measurements (Uttley et al. 2014). The 
iron line is the strongest feature aside of the continuum, 
see Fig. and it can potentially provide quite detailed 
information on the spacetime geometry close to the BH 
candidate. This technique relies on fits of the whole re¬ 
flected spectrum, but the spin measurement (or possible 
tests of the Kerr metric) is mainly determined by the iron 
line. For this reason the technique is often called the iron 
line method. 


The shape of the line is primarily determined by the 
background metric, the geometry of the emitting region, 
the disk emissivity, and the disk’s inclination angle with 
respect to the line of sight of the distant observer. In 
the Kerr background, the relativistic emission line pro¬ 
file emitted by an accretion disk illuminated by an X-ray 
corona with arbitrary geometry is typically parametrized 
by the BH spin a*, the inner and the outer edges of the 
emission region Tin and Tout, and the viewing angle i. The 
local spectrum Ig (the power radiated per unit area of 
emitting surface per unit solid angle per unit frequency) 
is determined by the reflection processes and the geome¬ 
try of the system. Assuming axisymmetry, A depends on 
the photon energy, the emission radius, and possibly on 
the emission angle (the angle of propagation of the pho¬ 
ton with respect to the normal to the disk) if the emission 
is not isotropic. The local spectrum may be modeled as 
a power-law A oc r~'^, where the emissivity index g is a 
constant to be determined by the fit. A more sophisti¬ 
cated choice is to assume an intensity profile A cx: for 

r < Tbreak and Ig oc r~^ for r > rbreak: where g = 3 corre¬ 
spond to the Newtonian limit in the lamppost geometry 
at large radii far from the X-ray source. With this choice, 


we have two free parameters for the emissivity profile, g 
and Tbreak- If we assume a point-like hot corona located 
on the axis of the accretion disk and just above the BH, 
it is possible to compute the emissivity index g = g(r) 
according to the height of the corona h (Dauser et al.| 


2013). Galactic BH binaries have line of sight velocities 


of just a few hundred km/s, which makes their relative 
motion negligible in the spectrum. In the case of AGN, 
the cosmological redshift of the source can instead be 
important for some objects and it can be an additional 
parameter of the model. 

While the continuum-fitting method requires indepen¬ 
dent measurements of the mass, the distance, and the 
viewing angle of the source, the iron line analysis does 
not require them. Mass and distance have no impact on 
the line profile, as the physics is essentially independent 
of the size of the system, while the viewing angle can be 
inferred during the fitting procedure from the effect of 
the Doppler blueshift. 

As in the continuum-fitting method, the iron line mea¬ 
surement assumes that the accretion disk can be de¬ 
scribed by the Novikov-Thorne model, and crucially de¬ 
pends on the inner edge of the disk being located at the 
ISGO radius. The latter point can be tricky for the iron 
line analysis, because there is controversy over this as¬ 
sumption in hard states, where “hard” refers to the fact 
that the hard power-law component is the dominant fea¬ 
ture in the spectrum. These states are those most fre¬ 
quently employed in reflection modeling. Specifically, it 
is argued whether the disk is in fact truncated at some 
radius larger than the ISGO (Parker et al. [2015 Plant 


et ah, 2015). At present, the evidence indicates that if 


truncation is present, it is likely to be mild (a factor of a 
few). 

In the case of stellar-mass BH candidates, one usually 
selects sources with high luminosity or sources in the 
hard-intermediate state, in which there are indications 
that the inner edge of the disk is closer to the compact 


object than in the hard states at low luminosities (Plant 
et al.[[20T^. In the case of supermassive BH candidates, 


this technique is more widely used, see the third column 
in Tab.|ll] In part, this is because we never have a precise 
estimate of the Eddington scaled accretion luminosity, 
due to the large uncertainties in the bolometric correc¬ 
tion and mass estimates. This works both ways so that 
it is hard to say whether we are in the thin disk range 
or not. In part, it is because we cannot choose a differ¬ 
ent spectral state due to the much longer timescales of 
supermassive BHs than stellar-mass BHs. 

If one were to fit data for a system in which the disk 
were truncated, and make the usual assumption that the 
inner-edge was at the ISGO, the fit would then incorrectly 
underestimate the value of spin (relative to the resulting 
Kerr prediction if it were truncated at the ISGO radius). 
The same question is at play for using reflection to test 
the Kerr metric with data from faint hard-states. In the 
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FIG. 10 Impact of the model parameters on the iron line profile: viewing angle i (top left panel), emissivity index q (top right 
panel), spin parameter a* (bottom left panel), and JP deformation parameter £3 (bottom right panel). When not shown, the 
values of the parameters are: i = 45°, q = 3, a, = 0.7, and £3 = 0. The outer radius is rout = Hsco + 100 M. Eohs is the 
photon energy measured by the observer far from the BH. The vertical axis is the number flux of photons in arbitrary units. 
The iron line profiles are computed with the code described in Bambi (2013c I. 


case of too high accretion luminosities, the inner part 
of the disk may instead be geometrically thick, and this 
would lead to overestimate the BH spin (or, otherwise, 
get a wrong constraint on the deformation parameter). 
Some very high values of the BH spins reported in Tab.|TT] 
are thus to be taken with caution (Reynolds 2014). 

Moreover, this is not the end of the story. As clearly 
shown in Dauser et al. (2013) within the lamppost set¬ 
up, a spin measurement is only possible when the corona 
is close to the BH candidate. In particular, the charac¬ 
teristic low energy tail of the iron line used to measure 
the inner edge of the disk can be produced only in the 
case of a compact corona close to a fast-rotating BH can¬ 
didate (Dauser et al. 2013). In the absence of the low 
energy tail, we may have either an extended corona far 
from the compact object or a slow-rotating BH candi¬ 
date. 


Bearing in mind these tricky points, the iron Ka line 
is used to measure the BH spin under the assumption 
of the Kerr background and it can potentially be used 
to test the Kerr metric. For a review, see e.g. [Reynold^ 


(2014). This technique can be used for both stellar-mass 
and supermassive BH candidates, because the iron line 
profile does not depend on the mass M. So it is currently 
the only available approach to probe the metric around 
supermassive BH candidates with the existing data. Cur¬ 
rent spin measurements of stellar-mass BH candidates 
with the iron line and under the assumption of the Kerr 
metric are reported in Tab. |lj A summary of spin mea¬ 
surements of supermassive BH candidates is shown in 

Tab. mg 

The impact of the model parameters on the iron line 
profile is shown in Fig. 10 where the emissivity profile 
has been modeled with a simple power-law le oc 


In some cases, there exist a few independent measurements of 
the spin of the same source in the literature. Measurements 
from different groups are often consistent, but sometimes they 
are not; the discrepancy is more likely due to systematic effects 
not fully under control (see the references in the last column for 
the details). 
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FIG. 11 contours with Nune = 10^ (top panels) and Nu^e = 10^ (bottom panels) from simulations of measurements of 
iron line profiles. The background geometry is the CPR metric with 63 and £3 as the only possible non-vanishing deformation 
parameters. The reference model is a Kerr BH with spin parameter a(. = 0.95 and inclination angle i' = 70°. In the left panels, 
we allow for a non-vanishing £3 while we assume £3 = 0. In the right panels we illustrate the converse measurement, namely 
£3 = 0 while £3 can vary. The ratio between the continuum and the iron line photon flux as well as the photon index of the 
continuum are also free parameters in the ht. The labels along the contour levels refer to the value of Ax^. See the text for 
more details. From Jiang et al. (20161. 


and q is constant. The iron line profile has a more 
complicated structure than the thermal spectrum of thin 
disks. However, as in the disk’s thermal spectrum, the 
main effect of the spin parameter and of possible devia¬ 
tions from the Kerr geometry is often similar: both a* 
and the deformation parameter change the ISCO radius. 
This in turn affects the extension of the low energy tail 
of the line. There is thus a parameter degeneracy, as 


shown in Johannsen & Psaltis (2013). As in the case of 


the continuum-fitting method, the technique can also be 


used to test the Kerr metric (Bambi 2013c Johannsen 


& Psaltis 2013 Lu & Torres 2003). 


The iron line is potentially a more powerful tool than 
the continuum-fitting method to test the Kerr metric, 
in the sense that in the presence of high quality data 
it is typically possible to get independent estimates of 
the spin and the deformation parameter, while with the 
continuum-fitting method high quality data may not be 
enough to break the parameter degeneracy. [Johannsen] 
& Psaltis (2013) estimated the required precision that 


observations with future X-ray missions have to achieve 
in order to measure potential deviations from the Kerr 
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metric with the iron line. They found that a precision of 
about 5% can constrain the absolute value of the JP de¬ 
formation parameter ea to be smaller than 1 if the source 
is a Kerr BH with spin parameter a* > 0.5 and the view¬ 
ing angle is 30° or 60°. The constraining power increases 
for fast-rotating Kerr BHs, because the inner edge of the 
disk is closer to the compact object. 

Jiang et al. (2015a|b 20161 presented a study on the 


possibility of constraining the JP and CPR deformation 
parameters. In the case of a bright AGN, a good ob¬ 
servation can have Mine ~ 10^ photons in the iron Ka 
line. In the case of a bright stellar-mass BH candidate, 
Mine may be up to two orders of magnitude higher, say 
Mine ~ 10®. However, in stellar-mass BH candidates the 
low energy tail of the iron line overlaps with the ther¬ 
mal component of the accretion disk, the ionization is 
higher and less easily modeled, and Compton broaden¬ 
ing plays a non-negligible role. All these effects make the 
brighter signal of stellar-mass BH lines messier to model. 
In Jiang et al. ( 2015a|b 20161, we simulated iron lines 
in Kerr and non-Kerr backgrounds to estimate plausi¬ 
ble constraints on the JP and CPR deformation parame¬ 
ters to be obtained from observations. We employed the 
standard analysis technique common throughout X-ray 
astronomy: the simulated spectra include Poisson noise 
and forward-folded fit after rebinninj^ to ensure that 
the count distribution in each spectral channel is well 
approximated by a Gaussian distribution, as required to 
use the as goodness-of-fit test. 

shows the Ay^ contours with Mine = 10® (top 


Fig. 


11 


panels) and Mine = 10"^ (bottom panels) assuming as a 
reference model a Kerr BH with spin parameter a* = 
0.95 and observed with a viewing angle i = 70°. The 


Forward-folded fitting is the term commonly used in X-ray as- 
tronomy to generically indicate the fitting procedure. The ac¬ 
tual spectrum measured by an instrument (in units of counts 
per spectral bin) can be written as 

C(h) = T J R{h,E)A(E)s{E)dE , (24) 

where h is the spectral channel, r is the exposure time, R(h, E) 
is the redistribution matrix (essentially the response of the in¬ 
strument), A{E) is the effective area, and s{E) is the intrinsic 
spectrum of the source. In general, the redistribution matrix 
cannot be inverted, and therefore it is not possible to obtain the 
intrinsic spectrum s{E). Forward-folded fitting means that one 
assumes a theoretical model for s{E) with a set of parameter 
values, convolves the model with the instrument response, and 
its result (called the folded spectrum) is compared with the ob¬ 
served spectrum through a goodness-of-fit statistical test. The 
process is repeated by changing the parameter values to find the 
best fit. For example, this approach is not necessary at optical 
wavelengths, where the redistribution matrix can be inverted, 
and therefore one can directly fit the theoretical model with the 
observed intrinsic spectrum of the source. Rebinning means that 
one groups data with a low count number to increase the num¬ 
ber of counts per bin. For more details, see, e.g.,|Arnaud et al.l 


spacetime geometry is described by the CPR metric and 
we have five free parameters, namely the spin a*, the 
deformation parameter (e^ or € 3 ), the viewing angle i, 
the ratio between the continuum and the photon iron line 
flux, say K, and the photon index of the continuum T. 
y^ is thus minimized over z, AT, and T. In the left panels, 
£3 is a free parameter and £ 3 = 6 . In the right panels, 
£3 = 0 and £3 can vary. It is evident that £| is relatively 
easy to constrain, while £3 is much more difficult. Even 
in the case of a fast-rotating Kerr BH observed with a 
large inclination angle, it is impossible to constrain £ 3 , in 
the sense that large positive values of £3 cannot be ruled 
out. The difference in the constraints on £3 and £3 can be 
understood noting that e\ strongly affects the position of 
the ISCO. The impact of £3 on the ISCO is weaker (it 
enters only in gtif,)- £3 mainly affects the propagation of 
the photons in the spacetime, but the impact on the time- 
integrated iron line measurement is weak. More details 
on the constraints on the CPR deformation parameters 


from the iron line can be found in Jiang et al. (2015b 


2016). 


As the photon count number increases, the constraints 
become stronger. It is not clear what the best current 
observations can provide and work on real data is in 
progress. It is however important to note that current 
theoretical models are quite simple. The emissivity pro¬ 
file, for instance, is often modeled with two power-law 
indices and a breaking radius, which is clearly an approx¬ 
imation. Even if we presently had excellent data with a 
large number of photons in the iron line, it would be 
presumably impossible to get reliable constraints on the 
Kerr metric. More realistic modeling would be necessary 
to prevent systematic effects from becoming dominant. 
In the end, any measurement is as good as its theoretical 
model, and current iron line models are phenomenologi¬ 
cal and oversimplify the astrophysical picture. 


1. Comparison between continuum-fitting and iron line 
measurements 


In the case of stellar-mass BH candidates, there are 
some objects in which the spin has been measured by 
both the continuum-fitting and the iron line methods, 
and the estimate is usually consistent, see Tab. |T] Let 
us note, however, that this is not enough to claim that 
these objects are Kerr BHs, because both techniques are 
mainly sensitive to the position of the inner edge of the 
disk and even in the presence of a non-Kerr background 
they should provide a very similar estimate of the spin pa¬ 


rameter when it is assumed to be the Kerr metric (Bambi 


2013d). The agreement between the two techniques is 


however a good result to believe in the robustness of 
the two measurements. The disagreement in the case of 
some objects is surely due to systematic effects and can¬ 
not be solved by postulating a non-Kerr metric. Only in 
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the case of very precise measurements, depending on the 
actual metric around the compact object, can the dif¬ 
ference between the continuum-fitting and the iron line 
measurements be attributed to deviations from the Kerr 


an extended low energy tail and the iron line of non¬ 
rotating wormholes looks like that from a fast-rotating 


Kerr BH (Bambi 2013e). An interesting feature of these 


iron lines is the peak at low energies. Like in the exotic 


solution (Bambi 2013d). 


compact stars in Bambi & Malafarina (2013), it is due to 


2. Exotic candidates 


the photons emitted at very small radii, which can more 
easily escape to infinity than in the case of a BH space- 
time. Can this feature be used as a smoking gun to rule 



out these spacetimes? It depends on the quality of the X- 

Note that the iron line can immediately rule out some 

ray data and on the specific metric (Zhou et al. 

2016 

). In 

exotic BH alternatives, or at least constrain their prop- 

Fig.[^ the iron line profile has been computec 

assuming 


erties. This is possible because in some spacetimes the 
iron line profile should have qualitatively different fea¬ 
tures that are not observed in real data or, vice versa, 
their iron line cannot mimic that observed in real data. 
This permits one to conclude that these objects cannot 
be the BH candidates of the Universe. 

Exotic compact stars represent a large class of BH al¬ 
ternatives. One usually assumes that general relativity 
holds, but the mass of the compact object can exceed the 
famous 3 Mq limit of a neutron star thanks to special 
equations of state. Boson stars also belong to this group 


of objects (Colpi et al. 1986 Jetzer, 1992) and they may 
also have a mass of millions or billions Mq , so that they 
may be the supermassive BH candidates at the center 
of galaxies. These spacetimes have typically no ISCO, 
namely equatorial circular orbits are always stable. One 
can set the inner edge of the disk at some radius r-m. 
Since the gravitational redshift of the photons emitted in 
the inner part of the accretion disk is never very strong, 
the iron line profile cannot have the characteristic low 
energy tail of the iron line expected in the case of very 
fast-rotating BHs and observed in the spectra of several 


BH candidates (Bambi & Malafarina 2013). In the case 


of fast-rotating BHs, the radiation emitted from the inner 
part of the accretion disk is strongly redshifted. This pro¬ 
duces an extended low energy tail in the iron line. In the 
case of compact stars, the gravitational redshift is usu¬ 
ally not very strong even at very small radii. Examples of 
iron line profiles from the accretion disk of exotic compact 
stars are shown in the left panel in Fig.[^ See [Bambi &:| 
Malafarina (2013) for more details. This argument does 


not imply the existence of a horizon, and consequently 
does not completely rule out scenarios with exotic com¬ 
pact stars, but at least it requires that the gravitational 
redshift close to BH candidates is extremely strong. It 
is not far from that expected in the vicinity of an ap¬ 
parent/event horizon, and this is not the case in typical 
exotic star models. 

Another example of strange candidates is repre¬ 
sented by the family of traversable wormholes discussed 


in Bambi (2013e). Their iron line profiles are shown in 
the right panel in Fig. [T^ These objects have no event 
horizon, like the compact stars mentioned in the previous 
paragraph. Nevertheless, their iron line profile has a low 
energy tail. Here it is actually impossible not to have 


an emissivity ~ corresponding to the Newtonian 
limit at large radii in the lamppost geometry. For such an 
emissivity profile, the iron lines in the Kerr background 
have no similar peaks, see Fig. jlOj However, in a cor¬ 
rect relativistic lamppost geometry one should expect a 
much steeper emissivity function at smaller radii and for 
a source close to the compact object, as a result of strong 
light bending (Dauser et al. 2013). Observations seem 


to require a high, or even very high, value of q at small 
radii (Wilkins & Fabian, 2011,2012 Fig. 13 shows the 


iron line in a Kerr spacetime with a* = 0.95 and i = 45° 
where the emissivity function is 


j ifr<rbreak, 

1 if ?■ > ’’break , 


(25) 


with Tbreak = 5 M and q assumes different values. As 
shown in Fig. |13[ such an emissivity profile could repro¬ 
duce a peak at lower energies! The reason is that the 
emissivity is much higher at small radii than the simple 
case r~^. Even if in the Kerr metric it is more difficult to 
escape from small radii, a much higher emissivity at small 
radii can balance the BH photon capture and have an iron 
line that looks like that of a traversable wormhole with 
lower emissivity at small radii. Eventually, only in the 
presence of high quality data is it possible to distinguish 
the difference between the astrophysical model and the 


metric and therefore between wormholes and BHs (Zhou 


et al., 2016). 


C. Iron line reverberation 

Within the corona-disk model with lamppost geome¬ 
try, the activation of a new flaring region in the corona 
illuminates the accretion disk and generates a time- 
dependent iron line profile due to the different propa¬ 


gation time for different photon paths (Reynolds et al. 


1999). Reverberation refers to the iron line signal as a 


function of time in response to a flash of radiation from 
the corona. The resulting line spectrum as a function 


However, this interpretation was criticized in jSvoboda et al.| 
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Eobs (keV) 
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FIG. 12 Iron line profiles produced from the accretion disks around a class of exotic compact stars (left panel) and a certain 
family of traversable wormholes (right panel). Eohs is the photon energy measured by the observer far from the BH. The 
vertical axis is the number flux of photon s in arbitrary units. See the text for more details. The iron line p rofil es are computed 
with the code described in Bambi (2013c I, employing the metrics discussed in Bambi & Malafarina ( 2013[) and Bambi (2013el. 



FIG. 13 Examples of iron line profile from a Kerr BH with 
spin parameter o* = 0.95 and observed from an inclination 
angle i = 45°. The emissivity is h oc (rbreak/r)® for r < 
rbreak = 5 M and h oc (rbreak/r)® for r > rbreak. See the text 
for more details. The iron line profiles are computed with the 
code described in Bambi (2013c I. 


of both time and across photon energy is called the 2D 
transfer function, which is related to fundamental prop¬ 
erties of the BH candidate and of the system geometry. 

In the case of supermassive BH candidates, iron line 
reverberation measurements are currently possible and 
probably represent the most convincing argument in sup¬ 
port of the fact the iron line originates from the inner 
part of the accretion disk, and therefore that its shape 
is determined by relativistic effects. For a review, see 
e.g. Uttley et al. (2014). However, with current X-ray 
facilities, because of the limited count rates in the iron 
line, reverberation measurements are not more powerful 
than the standard time-integrated ones in probing the 
spacetime around supermassive BH candidates. Future 


detectors with larger effective areas should be able to 
better study the temporal change in response to the ac¬ 
tivation of new flares and reverberation measurements 
can probably become a more powerful technique than 
time-integrated observations. 


Iron line reverberation mapping in a non-Kerr back¬ 
ground was investigated in Jiang et al. (2015a 2016). 
Fig. shows the constraints on the CPR deformation 
parameter Cg and eg from simulations. These plots can 
be directly compared with those in Fig. 11 for a time- 
integrated iron line measurement. The reference model 
is the same, namely a Kerr BH with spin parameter 
a* = 0.95 and a viewing angle i = 70°. The top pan¬ 
els refer to observations with a photon number count in 
the iron line Mine = 10^, and the bottom panels refer 
to the case Mine = 10"^- In the left panel, can vary 
and Cg = 0 is frozen. In the right panels, we have the 
opposite cases, so e| = 0 and Cg is free. The height of 
the source, the ratio between the continuum and the iron 
line photon flux, and the photon index of the continuum 
are also left as fit parameters. 


The case Mine = 10^ roughly corresponds to a good 
observation of a bright AGN with the current X-ray facil¬ 
ities. The reverberation measurement can constrain the 
background metric better than the time-integrated ob¬ 
servation. However, it is still problematic to measure the 
CPR deformation parameter eg. Reverberation mapping 
becomes much more powerful than the measurement of 
the time-integrated iron line when Mine increases. This 
is perfectly understandable. With a low photon number 
count, the large number of channels in the reverberation 
approach dilutes the photon count per channel and the 
intrinsic noise of the source frustrates the additional time 
information in the measurement. In the simulations with 
Mine = 10^, we find that the reverberation measurement 
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FIG. 14 Compare to Fig. 11 Ax^ contours with Nune = 10^ (top panels) and A^une = 10"* (bottom panels) from the analysis 


of the 2D transfer function. The reference model is a Kerr BH with spin parameter a( = 0.95 and inclination angle i' = 70°. 
In the left panels, we allow for a non-vanishing and we assume £3 = 0. In the right panels we consider the converse case, 
namely £3 = 0 and £3 can vary. The height of the source, the ratio between the continuum and the iron line photon flux, and 
the photon index of the continuum are also left as fit parameters. The labels along the contour levels refer to the value of Ax^. 
See the text for more details. FromljTang et aL](|2016[). 


can constrain the CPR deformation parameter £ 3 . In¬ 


terestingly, as discussed in Jiang et al. (2016), this is 


true even if the source is a slow-rotating Kerr BH ob¬ 
served from a low inclination angle. This is probably 
possible because £3 affects the photon propagation in the 
background metric. Time-integrated observations are not 
very sensitive to it, while the time information of the pho¬ 
ton propagation allows one to constrain £ 3 . 

Time-resolved measurements may also be possible from 


the observation of an X-ray AGN eclipse (Risaliti et al. 


2011). In principle, one could still exploit the variability 


of the source in order to probe different regions of the 


spacetime at different times, and better separate the rel¬ 
ativistic effects occurring near the BH candidate. How¬ 
ever, Cardenas-Avendano et al. (2016a) showed that this 


is not the case and an eclipse measurement does not have 
the advantage of a reverberation one. The difference be¬ 
tween the two approaches is related to the capability of 
separating photons from different parts of the disk. In 
the reverberation approach, photons emitted from differ¬ 
ent regions are detected at different times. In the eclipse 
scenario we have the opposite case, namely one studies 
the properties of the radiation from every region of the 
accretion disk from the non-detection of the photons from 
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that patch. 


D. Quasi-periodic oscillations 


Quasi-periodic oscillations (QPOs) are a common fea¬ 
ture in the X-ray flux of stellar-mass BH candidates. 
They appear as peaks in the X-ray power density spectra 
of the source. They are thought to be a very promis¬ 
ing tool for the future to get precise information on the 
spacetime geometry around BH candidates. However, 
there is currently no consensus on which mechanism is 
responsible for their production or even if it is a sin¬ 
gle mechanism or multiple ones. Many scenarios have 
been proposed, and in particular there are relativistic 


precession models (Motta et al. 2014 Stella & Vietri 


19991, diskoseismology models (Perez et al. 19971, res¬ 


et al. 2005), p-mode oscillations of accretion tori (Rez 


onance models (Abramowicz & Kluzniak 2001 Tbrok 


zolla et al., 2003). Other mechanisms may also be pos¬ 


sible. Interestingly, in most scenarios the frequencies of 
the QPOs are directly related to the characteristic or¬ 
bital frequencies of a test-particle (orbital or Keplerian 
frequency radial epicyclic frequency Vr, and vertical 
epicyclic frequency i^g), which are only determined by 
the background metric and are thus independent of the 
complicated astrophysical processes of the accretion. For 
this reason, QPOs are a promising technique to probe 
the metric around BH candidates. Moreover, the fre¬ 
quencies of QPOs can be measured with high accuracy, 
and therefore they can potentially provide more precise 
measurements than other techniques like the continuum¬ 
fitting and the iron line methods. 

Different models provide a different measurement of 
the parameters of the background metric, which means 
that QPO data cannot be used to test fundamental 
physics at this time. However, there is already some 
work exploring the possibility of using QPOs to test the 


Kerr metric (Aliev et al. 

2013 Bambi| 

2012b 2015a 

Jo- 

hannsen & Psaltis| 2011b 

Maselli et al. 

2015||Stuchlf 

c L 


Kotrlova, 2009). Even if QPO data can potentially pro¬ 


vide very accurate measurements, there is a fundamental 
degeneracy among the spin parameter and possible devi¬ 
ations from the Kerr solution, so one can typically only 
obtain a narrow allowed region on the spin parameter - 
deformation parameter plane. 


Johannsen & Psaltis (2011b) discussed the constrain¬ 


ing power of the diskoseismology model and of the 1:2 
resonance model involving the Keplerian and the radial 
epicyclic frequencies in the framework of the non-Kerr 


metric of Glampedakis & Babak (2006). In the disko¬ 


seismology scenario, the pair of high-frequency QPOs 
observed in the X-ray flux of some BH binaries can be 
identified as the lowest order gravity (g-modes) and cor¬ 
rugation modes (c-modes). Within the Kerr metric, the 
measurements of the two frequencies would provide the 


values of the BH mass and spin because the lowest or¬ 
der modes would occur near the ISCO, so the radius is 
fixed. When we want to test the Kerr metric and we have 
a non-vanishing deformation parameter, we need an in¬ 
dependent measurement of the mass, and we can get a 
value for the spin and the deformation parameter. In the 
case of the 1:2 resonance model, there is a degeneracy 
among the mass, the spin, and possible deviations from 
Kerr. A possible non-vanishing deformation parameter 
can only be constrained in the presence of independent 
measurements of the mass and the spin of the compact 
object. 


In Bambi (2012b), I considered a number of differ¬ 
ent resonance models to constrain the JP deformation 
parameter £ 3 . Consistently with Johannsen & Psaltis 


( 2011 b), any resonance model can measure just one num¬ 
ber of the spacetime geometry. In Bambi (2012b), the 
constraints of the resonance model for the BH candidates 


GRO J1655-40, XTE J1550-564, and GRS 1915-kl05 
were combined with the dynamical measurement of their 
mass and with their constraints from the continuum¬ 
fitting method. Assuming the Kerr metric, no resonance 
model provides a spin measurement consistent with the 
continuum-fitting method for the three BH candidates 
at the same time. In the case of a non-vanishing £ 3 , the 
3:1 resonance model involving the two epicyclic frequen¬ 
cies can be consistent with the continuum-fitting method 
estimates for the three objects. 


In Bambi (2015a), I studied the constraints on the 
JP deformation parameter £3 for the BH candidate in 
GRO J1655-40 assuming the relativistic precession model 
of Motta et al. (2014). GRO J1655-40 is a special source 
because it is the only object for which we have a de¬ 
tection of three simultaneous QPOs. Assuming that the 
three QPOs are associated with oscillations of the fluid 
flow at the same radial coordinate, in the Kerr metric 
there are three unknown quantities (the BH mass, the 
BH spin, and the radial coordinate of the fluid oscilla¬ 
tion) and it is possible to solve the system. In the case of 
a non-Kerr metric with only one non-vanishing deforma¬ 
tion parameter, we need an independent measurement of 
the BH mass. In the case of GRO J1655-40, there are two 


main measurements (Beer & Podsiadlowski, 2002 Shafee 


et al. 2006), which are not consistent each other. If we 


adopt the mass estimate of Beer & Podsiadlowski (2002), 
the constraint on £3 is consistent with the Kerr metric, 
but the estimate of the BH spin is not consistent with 
that from the continuum-fitting method. If we choose 
the mass estimate in Shafee et al. (2006), the measure¬ 
ments of the continuum-fitting method and the QPOs 
may be consistent for a non-vanishing £3. 
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FIG. 15 Polarization degree (left panel) and polarization angle (right panel) as a function of the photon energy for a 
Schwarzschild BH (solid l i nes) a nd a Kerr BH with spin parameter ajM = 0.9 (dashed lines) and a viewing angle i = 45°, 60°, 
( 2015|), under the terms of the Creative Commons Attribution License. 


and 75°. From Liu et al. 


E. X-ray polarization 


solution. 


The measurement of the polarization of the thermal 
radiation of thin accretion disks may become a new tech¬ 
nique to study stellar-mass BH candidates. This kind 
of measurement is not possible today, because there are 
no X-ray polarimetric missions. However, it will hope¬ 
fully be possible in the near future, for instance with the 
missions eXTP, XIPE, IXPE, and PRAXYS. 

The thermal radiation of a thin accretion disk is ini¬ 
tially unpolarized, but it gets polarized at the level of a 
few percent due to Thomson scattering of X-ray radia¬ 
tion off free electrons in the disk’s atmosphere. Because 
of relativistic effects (light bending and non-trivial par¬ 
allel transport in curved spacetime) more pronounced in 
the vicinity of the compact object, the degree and the an¬ 
gle of polarization of photons generated in the inner part 
of the accretion disk deviate from the Newtonian predic¬ 
tions, see Fig. |15| Assuming the Kerr background. X-ray 
spectropolarimetric observations of the thermal compo¬ 
nent could provide a measurement of the BH spin and of 
the inclination angle of the disk with respect to the line 


of sight of the observer (Li et al. 2009 Schnittman & 


Krolik 2009). 


As in the previous techniques, even the polariza¬ 
tion measurement may be used to test the Kerr met¬ 
ric (Krawczynski 2012 Liu et al. 2015). At present, 
there are only some preliminary studies about its con¬ 
straining capabilities. In Liu et al. (2015), we found that 
a polarization measurement cannot test the Kerr met¬ 
ric better than the continuum-fitting method, and it is 
definitively worse than high quality data of the iron line. 
The problem is still the strong correlation between the es¬ 
timate of the spin and possible deviations from the Kerr 


IV. THE SPECIAL CASE OF SgrA* 


As discussed in the previous section, the main prob¬ 
lem to test the Kerr metric is the parameter degeneracy. 
The spectrum of a Kerr BH can be usually reproduced 
quite well by non-Kerr objects with different values of 
the model parameters. To break the parameter degener¬ 
acy, it is usually helpful to have different measurements 
of the same BH candidate. If these measurements are 
sensitive to different relativistic effects, we may combine 
the observations and constrain possible deviations from 
the Kerr geometry. 

SgrA*, the supermassive BH candidate at the center 
of the Galaxy, may soon become quite an ideal object to 
test the Kerr metric. While there are currently no obser¬ 
vations suitable to test this BH candidate, we expect a 
number of unprecedented data with new facilities in the 


near future {Falcke & Markoff 2013). The combination of 


these measurements is a very promising approach to test 


the nature of SgrA*, see e.g. Bambi (2015b); Johannsen 


(2012). A recent review on tests of the Kerr metric with 


SgrA* is Johannsen (2015). However, at present we do 


not know if SgrA* has all the features to be an optimal 
source for testing the Kerr metric (e.g. high spin param¬ 
eter and large viewing angle). 


A. Black hole shadow 

The direct image of the accretion flow around a BH 
usually shows a dark area over a bright background. Such 
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FIG. 16 Ray-tracing calculations of the direct image of a Kerr 
BH surrounded by an optically thin emitting medium. The 
dark area is the BH shadow and its boundary corresponds 
to the apparent photon capture sphere. From [Falcke et al.| 
(20001. ©AAS. Reproduced with permission. 



FIG. 17 The function R{<j>) is defined as the distance between 
the center G and the boundary of the shadow at the angle 
0 as shown in this picture. See the text for more details. 
From Ghasemi-Nodehi et al. (20151, under the terms of the 
Creative Commons Attribution License. 


the accretion structure and on the emission mechanisms, 
the boundary of the shadow is only determined by the 
spacetime metric and the viewing angle of the observer. 
An accurate measurement of the direct image of the ac¬ 
cretion flow around a BH candidate can thus test the 
spacetime geometry around the compact object. 

Sub-millimeter very long baseline interferometry 
(VLBI) facilities should be able to resolve the shadow 
of SgrA* in the next few years (Doeleman et al. 20081. 


While it is currently not clear the effects of inevitable as- 
trophysical complications and the level of accuracy that 
can be reached in the measurement of the boundary of 
the image of the photon capture sphere in the observer’s 
sky, there has been much work to explore the possibility 
of testing the Kerr metric with the detection of a shadow 


and to calculate the shadows of non-Kerr BHs ( 

Abdujab- 

barov et al. 2013 Amarilla & Eiroa 2012 

2013 

lAmariila 

et al., 20101 Atamurotov et al. 2013a|b 

Bambi 2013f 

Bambi et al.[ 2012[ |Bambi & Freese |2009 

Bambi & 

Yoshidal |2010 |Ghasemi-Nodehi et al.[ 20151 , 

ohannsen 

& Psaltisl 20101 Li & Bambi 2014al |Schee & Stuchlfk 

2009 Tsukamoto et al. 2014 Wei & Liu 

2013 

I. 


In order to infer the spin and the deformation param¬ 
eters from a possible precise detection of the shadow of 
a BH, it is necessary to have a formalism to describe the 


posed in 

dhasemi-Nodehi et 

al. 

(2015 

I and illustrated 

in 

Fig. 17 

a more sophisticated method was presented 

in 

Abdujabbarov et al. 

(2015 

)]. First, we hnd the “cen- 


ter” C of the shadow. Its Cartesian coordinates on the 
image plane of the observer are 


J p{X,Y)XdXdY 
^ ~ "Yp{X,Y)dXdY^ 
y Jp{X,Y)YdXdY 
^ J p{X, Y)dXdY 


(26) 


where p{X, Y) = 1 inside the boundary of the shadow 
(which is a closed curve) and p{X, F) = 0 outside. As¬ 
suming a reflection-symmetric spacetime, the shadow is 
symmetric with respect to the X-axis and we can de¬ 
fine R{0) as the shorter segment between C and the 
shadow boundary along the X-axis. Defining the angle 


a dark area is commonly called the BH shadow (Falcke 


et al. 20001, even if the name may be a little bit mislead¬ 


ing. If the BH is surrounded by an optically thin emit¬ 
ting medium, the boundary of the shadow corresponds 
to the photon capture sphere as seen by a distant ob¬ 
server. The ray-tracing calculations of the direct image 
of a Kerr BH surrounded by an optically thin emitting 
medium are shown in Fig |16[ here the spin parameter is 
a* = 0.998 and the viewing angle is i = 45°. While the 
intensity map of the image depends on the properties of 


In what follows, strictly speaking we consider the boundary of 
the image of the photon capture sphere in the observer’s sky, 
which has a well defined border. If one fires a photon inside such 
a boundary, the photon is captured by the BH; if the photon is 
fired outside, it can then escape to infinity. In order to use such a 
description to analyze real data, it would be necessary to define a 
precise link between the boundary of the dark area in the direct 
image (real or simulated) of the accretion flow and the boundary 
of the image of the photon capture sphere in the observer’s sky. 
The intensity map of the image depends on the properties of the 
accretion flow and of the observational facilities (in particular on 
the wavelength of the observation), and there are no such data 
available today. 
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FIG. 18 BH shadows and R functions for CPR BHs with different values of the spin parameter a*, the deformations parameters 
£3 and £ 3 , and the inclination angle i. From Ghasemi-Nodehi et al. (20151, under the terms of the Greative Gommons Attribution 
License. 


(j) as shown in Fig. is the distance between the 

point C and the boundary at the angle 4>. The function 
R{(j))/R{Q) completely characterizes the shape of the BH 
shadow. Here we consider R{(j))/R{0) instead of R{4>) be¬ 
cause the latter cannot be measured with good precision, 
as it would require an accurate measurement of the dis¬ 
tance and of the mass of the BH, which is not available 
at the moment. Even the exact positions of the shadow 
on the image plane of the observer cannot be used to 
test the Kerr metric, because it is difficult to precisely 
identify the center X = Y = 0 of the source. 

Fig. [T8| shows some examples of shadows of CPR BHs 
and the associated R{(j))/R{0) function. In the top pan¬ 
els, the spin parameter is a* = 0.5, while in the bottom 
panels it is a* = 0.9. The inclination angle is always 
i = 85°, which is high and can thus maximize the rela¬ 
tivistic effects. In the left panels, Cg changes and £3 = 0 is 


frozen. In the right panels, we have the opposite case and 
£3 = 0 and £3 can vary. It is evident that mainly affects 
the size of the shadow, which increases (decreases) if £3 
decreases (increases). £3 alters the shape of the shadow 
on the side of corotating orbits, while there are no ap¬ 
preciable effects in the other parts of the boundary of 
the shadow. The peculiar boundary appearing for £3 = 2 
and 5 in the bottom right panel is due to the non-trivial 
horizons of these BHs (Bambi & Modesto 2011 1 . 


Note that VLBI observations do not directly image the 
accretion flow, and therefore they cannot directly mea¬ 
sure the shape of the boundary of the shadow discussed 
in this section. They instead sample the Fourier space 
conjugate to the sky image at a finite number of points. 
The boundary of the shadow can be obtained after im¬ 
age reconstruction. However, as mentioned, it is not clear 
whether a precise determination of the boundary of the 
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FIG. 19 Posterior probability density of spin and deformation 
parameters, marginalized over all other quantities, from the 
existing EHT data. The solid, dashed, and dotted white lines 
correspond, respectively, to the 1-, 2- , and 3-a boundaries. 
The dashed gray lines correspond (from top to bottom) to 
spacetimes with risco/M = 6, 5, and 4. The grayed region in 
the lower right corresponds to metric with risco/M < 4 and 
it is neglected because the calculations may be affected by 
some pathological properties of this non-Kerr metric. Here 
e is the deformation parameter of the non-Kerr metric of 


Glampedakis & Babak (20061. From Broderick et al. (20141. 
©AAS. Reproduced with permission. 


shadow at the level necessary to test the Kerr metric is 
eventually possible, because systematic effects may pre¬ 
vent it. 

The Event Horizon Telescope (EHT)p^is a project in¬ 
volving mm and sub-mm observatories equipped with 
VLBI instrumentation to get high resolution images of 
the accretion flow around supermassive BH candidates 
at 230 and 345 GHz (Fish et al. 20131. One of the 
main goals of this experiment is the observation of the 
shadow of SgrA*. The existing mm-VLBI observations 
have been done with only three stations (respectively in 
Hawaii, California, and Arizona). Employing a radia- 
tively inefficient accretion flow model, [Broderick et al.| 
(2014) have explored the capability of present observa¬ 
tions to constrain possible deviations from the Kerr ge¬ 


ometry. They used the quasi-Kerr metric of Glampedakis 


& Babak (2006). In their simulations, we cannot see the 


exact shape of the apparent photon capture sphere be¬ 
cause at the wavelengths accessible to mm-VLBI it is 
partially obscured by the optically thick structure on the 
approaching side of the accretion flow. The result of this 


analysis is shown in Fig. 19 


http://www.eventhorizontelescope.org/ 


a [mpc] 

0.1 1 



FIG. 20 Precession timescale from the mass (M), the spin 
(S), the quadrupole moment (Q), and stellar perturbation 
(P) for a pulsar orbiting SgrA* as a function of orbital period 
Pb, assuming an orbital eccentricity of 0.5 and 10^ objects, all 
with mass M = Mq, within 1 mpc around SgrA*. From[Llu| 


et al. (20121. ©AAS. Reproduced with permission. 


B. Accurate measurements in the weak field 


The nature of SgrA* can be potentially investigated 
even with accurate measurements of the spacetime met¬ 
ric at relatively large radii. In this case, the gravita¬ 
tional field is weak, M/r 1, and we can adopt an 
approach similar to the PPN formalism of Solar Sys¬ 
tem experiments. If a spacetime is stationary, axisym- 
metric, asymptotically flat, Ricci-flat outside the source, 
and analytic about the point at infinity, its metric in 
the region outside the source can be expanded in terms 


of mass moments and current moments Se (Geroch 
1970 Hansen 1974p In the case of reflection symme¬ 


try, the odd M-moments and the even S'-moments are 
identically zero, so that the leading order terms are the 
mass Mq = M, the spin angular momentum Si = J, and 
the mass quadrupole moment M 2 = Q. In the case of a 
Kerr BH, the metric is completely determined by M and 
J, and all the moments Mg and Se, are locked to the mass 


The expansion in multipole moments is also possible when the 
spacetime is not axisymmetric, but in this case the mass and 
the current moments of order i are tensors. If the spacetime is 
axisymmetric, there are some simplifications, and the mass and 
the current moments of any order £ are completely determined 
by two scalars, namely Mi and Si- In the case of tests of the 
Kerr metric in the weak field, it is common to assume that the 
spacetime is axisymmetric, because it sounds a physically plausi¬ 
ble hypothesis and simplifies the problem. Let us also note that 
some assumptions may not hold in some relevant cases. For in¬ 
stance, BH solutions in alternative theories of gravity may not be 
Ricci-flat, and an example is the case of BHs in Einstein-dilaton- 
Gauss-Bonnet gravity. The analyticity assumption may also not 
hold, as in the case of the presence of a massive scalar field with 
a Yukawa type solution. 
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and the spin by 

+ = , (27) 

where i is the unit imaginary number. In particular, the 
mass quadrupole term is Q = 

This approach has the advantage that it is quite gen¬ 
eral and relies on the small number of assumptions. Even 
the requirement of the Ricci-flat spacetime may approxi- 
matively hold in many cases. Here the spin measurement 
is really a spin measurement, and it should not be cor¬ 
related to possible deviations from the Kerr solution in 
the near horizon region because the latter corresponds 
to higher order corrections. Such a measurement could 
be combined with a measurement in the strong gravity 
regime, which is usually a constraint on the spin and pos¬ 
sible deviations from the Kerr geometry. In some favor¬ 
able cases, it could be possible to determine also the mass 
quadrupole moment Q. This could permit us to test the 
Kerr metric at the quadrupole term, because in the case 
of a Kerr BH one must recover Q = Higher or¬ 

der terms can unlikely be tested, but clean measurements 
of J and Q would instantly be very helpful if combined 
with observations in the strong gravity field. 


1. Radio pulsars 


It is thought that a population of radio pulsars is or¬ 
biting SgrA* with a short orbital period and there is al¬ 


ready an intense work to detect these objects (Lorimer 


& Kramer 20051. For instance, Chennamangalam & 


Lorimer (20I4| argue that there may be ~ 200 pulsars 


in the inner parsec region (orbital period ^ 10^ yrs). 
Accurate timing observations of a radio pulsar orbiting 
SgrA* in a very close orbit (.^ I yr) would allow a precise 
measurement of the mass, the spin, and - in exceptional 
cases - even of the mass quadrupole moment of the super- 
massive BH candidate at the center of our Galaxy if the 


system is sufficiently free of external perturbations (Liu 


et ah, 2012). 


Because of the high electron density in the ionized 
gas at the center of the Galaxy, this kind of observa¬ 
tion must be made at much higher frequencies than those 
normally used for pulsar timing, which further challenges 
these measurements. A more serious problem is the pos¬ 
sible presence of a population of stars or BHs orbiting 
very close to SgrA*. The presence of these bodies may 
strongly affect, or even prevent, a clean measurement of 
the spin and the quadrupole moment of SgrA* with the 
radio pulsar. At the moment it is impossible to make 
predictions about the potentialities of future observations 
because we do not know the actual situation near SgrA*. 

Assuming a population of 10^ objects with a mass 
M = Mq isotropically distributed within 1 mpc around 


SgrA*, Liu et al. (2012) estimated the necessary orbital 


period of the pulsar to get a measurement of the mass. 


the spin, and the quadrupole moment of SgrA*. Fig. 20 


shows the timescales of secular orbital precession for a 
pulsar orbiting SgrA* as a function of its orbital pe¬ 
riod: the contributions are from the mass monopole M 
(pericenter advance), the spin S (frame dragging), the 
quadrupole moment Q; and stellar perturbation P. The 
orbital eccentricity is assumed to be 0.5. The precession 
timescale of the pericenter advance is already lower than 
that of stellar perturbation for a 10 year orbital period, 
which means that the observation of a radio pulsar with 
an orbital period less than 10 years can already be used 
to estimate the BH mass M. The measurement of the 
spin would require an orbital period less than 0.5 years to 
have the contribution from frame dragging significantly 
above that from stellar perturbation. The measurement 
of the quadrupole moment Q requires an orbital period 
less than 0.1 years. These estimates have to be taken as 
a general guide and the actual situation may be different. 
For instance, a population of 10 Mq BHs in this region 
may completely spoil the measurement of the parame¬ 
ters of the metric around SgrA*, as well as a significant 
anisotropy in the distribution of these bodies may chal¬ 
lenge it. 

Assuming a not too optimistic situation, in which we 
can observe a radio pulsar in an orbit of several months, 
timing observations could measure the mass and the spin 
of SgrA*. Since the pulsar is in the weak field of SgrA*, 
this would be a clean measurement of the spin parame¬ 
ter a*, namely independent of the higher order multipole 
moments of the spacetime. First, such a measurement 
should satisfy the Kerr bound |a*| < 1, because other¬ 
wise SgrA* could not be a Kerr BH. Second, the spin 
measurement could be combined with other observations 
of the strong field (shadow, hot spot, etc.) in which there 
is typically a strong correlation between the estimate of 
the spin and possible deviations from the Kerr solution 
to break such a degeneracy. In this manner we could test 
the Kerr metric. 


2. Normal stars 


Even normal stars in compact orbits can be used to 
probe the weak gravitational field of SgrA* and measure 
its spin parameter and, possibly, its mass quadrupole mo¬ 


ment. The idea was proposed in Will (2008) and further 


explored in Angelil et al. (2010); Merritt et al. (2010). 


If SgrA* is rotating fast, the observation of at least 
two stars with an orbital period of O.I years or less and 
in orbits with a high eccentricity, say ^ 0.9, may pro¬ 
vide a measurement of the mass, the spin, and the mass 
quadrupole moment of SgrA* and thus test the Kerr na¬ 
ture of this object at the level of the quadrupole term. 
Today we know stars with an orbital period as short as 
10 years. These objects are still too far from SgrA* and 
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currently no relativistic effects in their orbit are observed 
(but they should be observed in the near future). How¬ 


ever, observational facilities like GRAVITY (Eisenhauer 


et al.l 20081, with the capability of accurate astrometric 


measurements at the level of 10 fj,as close to SgrA*, may 
observe stars with a sufficiently short orbital period to 
test the Kerr metric. 


Will (20081 proposed to test the Kerr metric by mea¬ 


suring the precession of the orbital plane of these stars. 
The advance of the pericenter of these stars is domi¬ 
nated by the mass term of the supermassive object, while 
the contribution of the spin and the quadrupole moment 
would be subdominant and difficult to estimate. On the 
contrary, in the weak field limit, the precession of the 
orbital plane is only determined by the spin (through 
the frame-dragging) and the mass quadrupole moment; 
see Will (2008) for the details. Here the dominant con¬ 


tribution comes from the spin, but in the case of suffi¬ 
ciently compact orbits it is also possible to infer the mass 
quadrupole moment Q. As in the pulsar case, a measure¬ 
ment of the spin could already be useful in combination 
with other measurements probing the strong gravity field. 
In the case a quadrupole moment measurement is also 
available, one can check a posteriori whether it satisfies 
the Kerr relation Q = 

Recently, Zhang et al.| ( 2015| ) showed that measure¬ 
ments of the spin of SgrA* are possible even with 
stars with orbital configurations similar to those already 
known, as in the case of long-term high precision obser¬ 
vations. 


C. Hot spots 


SgrA* exhibits powerful flares in the X-ray, near¬ 


infrared, and sub-millimeter bands (I 

)odds-Eden et al. 

2010 

Genzel et al. 

2003 

Trap et al. 

2011 

). During a 


flare, the flux increases up to a factor of 10. There are 
a few flares per day. Every flare lasts 1-3 hours and has 
a quasi-periodic substructure with a timescale of about 
20 minutes, see Fig. These flares may be associ¬ 
ated with blobs of plasma orbiting near the ISCO of 


the supermassive BH candidate (Hamaus et al. 2009), 


even if current observations cannot exclude other ex- 


planations ( 

Markoff et al. 

2001 

Tagger & Melia 

2006 

Yusef-Zadeh et al. 

2006 

). Temporary clumps of matter 


should indeed be common in the region near the ISCO 


radius (De Villiers et al. 2003) and, if so, they may be 


studied by the GRAVITY instrument for the ESO Very 


Large Telescope Interferometer (VLTI) (Eisenhauer et 
ahl [20081). 


The radiation emitted by a blob of plasma orbiting the 
strong gravity region of a BH candidate is significantly 
affected by the metric of the spacetime and can poten¬ 


tially be used to test the Kerr metric (Li et al. 2014). 


Fig. shows the spectrogram, namely the spectrum as 


a function of time, of a monochromatic blob of plasma 
orbiting the ISCO of a Schwarzschild BH and observed 
at a viewing angle of f = 60°. An accurate measurement 
of the spectrogram of a similar blob of plasma would 
be surely an ideal tool to test the metric around SgrA*. 
However, in the reality the situation is much more com¬ 
plicated. The astrophysical model (shape and size of the 
blob of plasma, spectrum of the blob of plasma in its 
rest-frame, etc.) is usually much more important than 
the small features associated with the relativistic effects 
characterizing the background metric (Li et al. 2014). 

At the moment, it is not clear if tests of the Kerr metric 
are possible with this approach. In some spacetimes, the 
photon capture radius can be significantly different from 
that of Kerr BHs, and in this case it is possible to identify 


specific signatures of these metrics (Li & Bambi 2014b 


Liu et al. 2015). In general, it seems more likely that 


relativistic effects cannot really be identified and eventu¬ 
ally the radiation from a blob of plasma can only provide 
an estimate of the orbital frequency. Since the observed 
period of the quasi-periodic substructure of the flares of 
SgrA* ranges from 13 to about 30 minutes, the orbital 
radius of these hot spots should change and be at a ra¬ 
dius larger than the ISGO. For a 4 million Mq Kerr BH, 
the ISGO period ranges from about 30 minutes (a* = 0) 
to 4 minutes (a* = 1 and corotating orbit). The shortest 
period ever measured is 13 ± 2 minutes, and it may be 
an upper bound for the ISGO period. In the Kerr met¬ 
ric, such a measurement translates into the spin estimate 
a* > 0.70 ± 0.11 (Trippe et al. 2007). 

In the case of a metric with a non-vanishing deforma¬ 
tion parameter, there is a degeneracy between the esti¬ 
mate of the spin and possible deviations from the Kerr 
solution. An example of the possible constraints is shown 
in Fig. |23[ Such a degeneracy may be broken with an¬ 
other measurement, for instance a precise estimate of the 
spin parameter by the observation of a radio pulsar. The 
latter would be independent of the deformation param¬ 
eter ea because the pulsar would probe the weak field 
limit, where the metric can be expanded in M/r. The 
leading order term in is subdominant with respect to 
the leading order term in a* = J/M'^ and the pulsar data 
may not be able to measure ea. 


D. Spectrum of the accretion structure 

Today, the two most popular and widely used tech¬ 
niques to probe the metric around BH candidates are 
the continuum-fitting and the iron line methods. Both 
approaches require geometrically thin and optically thick 
accretion disks and employ the Novikov-Thorne model. 
However, any accretion structure determined by the met¬ 
ric of the spacetime can potentially be used to test BH 
candidates. 

The accretion structure around SgrA* seems to be 
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FIG. 21 Light curve of a NIR flare of SgrA* with the charac¬ 
teristic quasi-periodic substructure with a timescale of about 
20 minutes. The arrows at the bottom indicate the peaks of 
the snbstrnctnre. Reprint with permission from |Genzel et J] 
( |2003[ ) and [Gillessen et al.| ( |2010[ ) . 







FIG. 22 Spectrogram of a hot spot orbiting a Schwarzschild 
BH at the ISCO radins from an inclination angle i = 60°. The 
vertical axis is the ratio between the photon energy measured 
by a distant observer and the photon energy at the emission 
point. The color scale indicates the energy flux (in arbitrary 
units). We see two tracks becanse one is the spectrogram of 
the primary image, the other one is for the secondary image. 
The points labeled 1-5 refer to another figure in the original 
paper and can be ignored here. From Li et al. (20141. 


a radiatively inefficient advection dominated accretion 


flow (Narayan et al. 

1995 

). There are a few differ- 

ent analytic models, with several variants, in the liter- 

ature ( 

Yuan & Narayan 

2014 

). In principle, one could 


use a model to fit the data and infer the free parameters 



FIG. 23 Assuming that the hot spot is at the ISGO radius, the 
measurement of its light curve can only select the spacetimes 
with the same orbital frequency. In this plot, the reference 
model employs a Kerr BH with spin parameter a* = 0.4. 
The allowed regions are those between the two blue dotted 
lines (Ay^ < 1 ) and the two red dashed lines (Ay^ < 10 ). 
The possible measurement of the spin through a binary pulsar 
would provide the constraint given by the yellow area, and the 
combination of the two observations could potentially break 
the degeneracy between a* and £ 3 . See the text for more 
details. From E eFaLldMIl ). 


of the model, including those related to the geometry 
of the spacetime. An explorative work in this direction 


was presented 

in 

Lin et al. 

2015 

1 , which employs the ion 

torus model of 

Straub et al 

(2012) 

Vincent et al. 

(2015). 


This approach currently faces a few problems. First, 
we do not know the correct model, and different mod¬ 
els presumably provide different results. Second, even 
simple analytical models have several free parameters 
that should be inferred by fitting the spectrum. The 
Novikov-Thorne model for thin disks follows from basic 
rules like the conservation of mass, energy, and angular 
momentum; the model is already quite constrained with 
a reasonable number of free parameters. The accretion 
structures necessary to describe the flow around SgrA* 
are more complicated and unconstrained. The advantage 
is that the spectrum of the accretion structure around 
SgrA* seems to have many features and this may break 
the degeneracy among the parameters of the model. For 
the moment, there is no accurate measurement of the 
full spectrum, and therefore it is impossible to constrain 
the model. However, this is presumably possible in the 
future. 
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V. OTHER APPROACHES 
A. Tests in weak fields 

SgrA* is not the only BH candidate for which we may 
get accurate measurements at relatively large radii. How¬ 
ever, in the case of other BH candidates we may not be 
so fortunate as to have independent measurements in the 
strong gravity field. In order to test the Kerr metric 
with these objects, we need to have measurements good 
enough to determine M, J, and Q. One can then check a 
posteriori whether Q = as expected in the case 

of a Kerr BH. 

The ideal candidate for this kind of test is a pulsar 
binary in which the companion is a stellar-mass BH can¬ 


didate (Wex & Kopeikin, 19991. A similar system is not 


known at the moment, but it should not be too rare and 
there is no reason to believe that it cannot be found in 
the future. It is also possible that the signal of a binary 
pulsar with a BH companion is already in the available 
radio data, but the data have not yet been analyzed. The 
identification of a new pulsar is a very time-consuming 
process. After the measurement of the period of the pul¬ 
sar, it is just an issue of time and accurate measurements 
of the system can be obtained thanks to the fact that a 
pulsar is like a precise clock. 

Even if the uncertainty is large in comparison with 
what could be possible with a pulsar binary, the mea¬ 
surement of the mass quadrupole moment of a BH can¬ 


didate has been done in Valtonen et al. (2010) and it 


could be relatively improved in the future (Valtonen et 


al. 2011). The object is the supermassive BH candidate 


in the quasar OJ287. Optical observations show a quasi- 
periodic light curve characterized by two timescales, one 
of ^12 years and another of ~60 years. The interpreta¬ 
tion is that the system is a binary BH, with the secondary 
BH orbiting the more massive primary one with an or¬ 
bital period of ~12 years and a periastron precession of 


60 years (Lehto & Valtonen 1996). The observed ma¬ 


jor outbursts occurring every ~12 years are thought to 
be due to the secondary BH that crosses the accretion 
disk of the primary. Within this interpretation, |Valtonen| 
et al. (2010) employed a 2.5 PN (Post-Newtonian) accu¬ 


rate orbital dynamics to fit current observational data. 
Since the mass quadrupole moment interaction term en¬ 
ters at the 2 PN order, it is possible to constrain the 
mass quadrupole moment of the primary BH. Writing 
the quadrupole moment as Q = where g = 1 

for a Kerr BH, current observational data provide the 


measurement (Valtonen et al. 2010) 


g = 1±0.3, 


(28) 


namely the mass quadrupole moment is tested at the 
level of 30%. In the next few years, this test could be 


improved at the level of 10% (Valtonen et al. 2011). 


B. Black hole jets 

Jets and outflows are quite a common feature of ac¬ 
creting compact objects. In the case of stellar-mass BH 


candidates, we observe two kinds of jets (Fender et al. 


2004). Steady jets are observed when the source is in 


the hard state. Transient jets typically appear when the 
source switches from the hard to the soft state. 

The actual mechanism responsible for the formation of 
these jets is currently unknown. One of the most popular 


scenarios is the Blandford-Znajek mechanism (Blandford 


& Znajek, 1977), in which magnetic fields threading the 


BH event horizon are twisted and can extract the rota¬ 
tional energy of spinning BHs producing an electromag¬ 
netic jet. Numerical simulations show that the process 
can be very efficient and strongly depends on the BH 


spin (McKinney 2005 Tchekhovskoy et al. 2010, 2011). 


The Blandford-Znajek mechanism is usually considered 
for the formation of steady jets. Other mechanisms do 
not involve the BH spin but still require magnetic fields 
to collimate the jets. 

Observations of a possible correlation between esti¬ 
mates of jet power and BH spin measurements are contro¬ 


versial. In Fender et al. (2010), the authors claimed that 


there is no evidence of a correlation between the jet power 
and the spin measurements of BH binaries reported in 
the literature with the continuum-fitting and the iron 


line methods. See also Russell et al. 


tails. Narayan & McClintock (2012) proposed that the 


(2013) for more de- 


Blandford-Znajek mechanism may be responsible for the 
formation of transient jets and showed that there is a 
correlation between jet power and the most recent spin 
measurements with the continuum-fitting method. The 
left panel in Fig. [24| shows the data used in [Fender et al. 
(2010) for steady jets and continuum-fitting spin mea¬ 


surements. The right panel shows the data reported in 


Narayan & McClintock 

(2012 

1 . 


The discrepancy between 

Fender et al. 

to 

o 

I—‘ 

o 

) and 

Narayan & McClintock 

(2012 

) can be easily understood. 


The two groups use different spin measurements and a 
different method to estimate the jet power. Moreover, 
as shown in Fig. |24[ we only have a few measurements 
with large error bars. In the future, with a larger number 
of measurements of the jet power and more precise spin 
estimates, it will be possible to test the existence of a 


correlation between jet power and BH spin, see Steiner 


et al. (2013). 


If the actual mechanism responsible for the formation 
of steady or transient jet is the Blandford-Znajek one, 
the measurement of the jet power could be used to esti¬ 
mate the BH spin if we assume the Kerr metric ([Steiner 


et al. 


2012 c 


2013) or to test the Kerr metric otherwise (Bambi 


d Pei et al. 2016). The key-point is that the 


estimate of the jet power is (typically) quite indepen¬ 
dent of the nature of BH candidates, while the measure¬ 
ment of the spin does depend on the background met- 
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ric. If the Blandford-Znajek mechanism is responsible for 
the formation of transient jets, as suggested in Narayan] 
& McClintock (20121, we can constrain possible devi¬ 
ations from the Kerr solution. If BH candidates were 
not Kerr BHs, the spin measurements would be differ¬ 
ent and we would lose the correlation between BH spin 
and jet power, see Bambi| (2012c). On the contrary, if 
the Blandford-Znajek mechanism is responsible for the 
formation of steady jets, the absence of a correlation 
between BH spin and jet power found in [Fender et al.| 
( 2010 ) can be explained with the fact that the estimate 


of the BH spin is wrong, because it was obtained assum¬ 


ing the Kerr metric. As shown in Bambi (2012d), if the 


continuum-fitting spin measurements are reanalyzed in a 
non-Kerr background it is possible to find a correlation 
between BH spin and jet power. At present these are 
just speculations based on a small number of data with 
large uncertainty, but they may provide some interesting 
results in the future. 


C. Evolution of the spin parameter 

As seen, there is a strong connection between the spin 
and possible deviations from the Kerr solution. In partic¬ 
ular, most observations are sensitive to both parameters 
and their measurement is usually correlated. In actuality, 
the connection between the spin and the deformation pa¬ 
rameters is even stronger. For instance, the Kerr metric 
with I a* I > 1 can unlikely have astrophysical implica¬ 
tions: as discussed in Subsection |H.C[ it is not clear if 
a similar object can be created by a physical mechanism 
and, even if it could be created, the spacetime would 
be very unstable and should quickly evolve to something 
else. The same considerations may hold for non-Kerr BHs 
and other exotic compact objects. Depending on the spe¬ 
cific spacetime metric and gravity theory, there may be a 
critical value of the spin parameter above which the sce¬ 
nario cannot be relevant in astrophysics, even in the case 
its gravity theory is right. Such a critical value may be 
either higher or lower than I, because |a*| = I does not 
correspond to any special property outside of the Kerr 
metric. 

Neglecting possible instabilities, which would depend 
on the unknown gravity theory, one may get a rough es¬ 
timate of the maximum BH spin as follows. Considering 
the astrophysical processes that can spin up and down 
a BH, the most efficient mechanism to create very fast- 
rotating bodies seems to be a prolonged accretion from 


a thin disk (Berti & Volonteri 2008). In the Novikov- 


Thorne model, the inner edge of the disk is at the ISCO 
radius. The gas’s particles reach the ISCO and then im¬ 
mediately plunge onto the compact object with specific 
energy Eisco and specific angular momentum TiscOj 
namely with their value at the ISCO radius. The mass 
M and the spin angular momentum J of the compact 


object change, respectively, by 

5M = EiscoSm, 
SJ = Tisco^w, 


(29) 

(30) 


where 6m is the gas rest-mass. The evolution of the 
spin parameter turns out to be governed by the following 
equation (Bardeen 1970) 


do* 


_ _ 1 Aisco 

d\iiM M Eisco 


2 o* 


(31) 


where here the small effect of the radiation emitted by the 
disk and captured by the object is neglected. Prolonged 
disk accretion is a very efficient mechanism to spin the 
compact object up till an equilibrium spin parameter oj'^. 


which is reached when the right-hand side of Eq. (31) 
becomes zero. If, for some reason, the spin is o* > oC, 
the accretion process spins the compact object down. 

In the case of the Kerr metric, it is possible to integrate 


Eq. (31) and find an analytic expression for the spin pa¬ 


rameter o* as a function of the BH mass M (Bardeen 


1970) 


a* = < 

[2 Mo 

V 3 M 

4- 

^18^-2 


U 




if M < VQMq , 
if M > V&Mo , 


(32) 


assuming an initially non-rotating BH with mass Mq. In 
Eq. (32), the equilibrium spin parameter is a*'^ = 1, 


which is reached after the BH has i ncrea sed its mass 
by a factor of ^/6 « 2.4 (Bardeen 1970). If we in¬ 


clude the effect of the radiation emitted by the disk 
and captured by the BH, we obtain the Thorne bound 


0.998 (Thorne 1974), because radiation with an¬ 


gular momentum opposite to the BH spin has a larger 
capture cross-section. 

One can repeat the calculations in a non-Kerr metric 
and find the corresponding equilibrium spin parameter 


which may be either larger or smaller than 1 (Bambi 
2011c|d ). In the case of the JP metric with non-vanishing 
£ 3 , the curve of a*'^ is the black solid line shown in Fig. 25 


If the object is on the left of the curve, an accretion disk 
spins it up to reach the equilibrium spin parameter. If the 
object is on the right of the curve, the accretion process 
spins the body down. As we can see from Fig. ^ > 1 

for £3 < 0 and a*"^ < 1 when £3 > 0. Accretion from thick 
disks might be a little bit more efficient, but the difference 


would be small (Li & Bambi 2013a). 


In the case of the supermassive BH candidates in galac¬ 
tic nuclei, the initial value of their spin parameter is ir¬ 
relevant, as these objects have increased their mass by 
several orders of magnitude and their spin parameter 
has evolved accordingly. Prolonged disk accretion is the 
most efficient mechanism to get a high spin and therefore 
we can conclude that these objects cannot have today 
a spin parameter larger than a*'^ ( Bambif 2011b). At 
the same time, it is possible to provide a lower bound 
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Steady Jets Transient Jets 




FIG. 24 Left panel: data reported in Fender et al. (20101 to claim the absence of evidence for a correlation between the 


power of steady jets and BH spin measurements of the continuum-fitting method. Right panel: data reported in|Narayan fc| 


McClintock ( 2012[) to claim the evidence for a correlation between the power of transient jets and BH spin measurements of 


the continuum-fitting method. 




FIG. 25 JP BHs with deformation parameter 63 . The black solid line corresponds to the equilibrium spin parameter for 
a thin disk, as it is inferred from Eq. (311. The red solid, green dashed, and blue dotted curves are, respectively, the contour 
levels of the Novikov-Thorne radiative efficiency = 0.15, 0.20, and 0.25. Supermassive BH candidates in galactic nuclei 
must be on the left of the black solid line and observations show that their radiative efficiency can be high or even very high. 
The allowed region is thus on the left of the black solid line and on the right of the r^NX contour levels. Th e right panel is the 
enlargement of the region of intersection of the curves with higher spin. See the text and the original paper Bambi (2011b I for 
more details. 


on their radiative efficiency, either as a mean radiative 


efficiency from the Soltan argument (Elvis et al. 2002 


Wang et al. 2006) or for specific sources (Davis & Laor 


2011). The radiative efficiency of these objects seems to 


be high, or even very high, but the exact estimate is more 
controversial. Fig. shows the contour levels of the 
Novikov-Thorne radiative efficiency = 0-15 (red solid 
curve), 0.20 (green dashed curve), and 0.25 (blue dotted 
curve). The observed high radiative efficiency constrains 
the spacetime geometry around supermassive BH candi¬ 


dates (Bambi, 2011e, 2012e|. Assuming the conservative 
bound ryNT > 0.15 and imposing that a* < aj'*, from 
Fig. we see that supermassive BH candidates are con¬ 


strained to be in the region bounded on the left by the red 
solid curve of ?7 nt = 0.15, and on the right by the black 
solid line of the equilibrium spin parameter. Moreover, 
it is possible to estimate an upper bound for the spin 
parameter of supermassive BH candidates, at the level 
of |a*| ^ 1.2 (1.1, 1.05) if one assumes the constraint 
??NT ^ 0.15 (respectively 0.20, 0.25), and this seems to 
be only very weakly dependent on the exact non-Kerr 
metric adopted for the calculations ([Bambi 2011bI. 


In the case of stellar-mass BH candidates, it is com¬ 
monly thought that the spin of the compact object is 


natal (King & Kolb 1999); however, also see Fragos & 


McClintock (2015). The point is that stellar-mass BH 
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FIG. 26 Evolution of the spin parameter a, (left panels) and of the Novikov-Thorne radiative efficiency ?7 nt = 1 — f^isco (right 
panels) as a function of the amount of matter accreted onto an initially non-rotating BH for four different initial BH masses. 
The top panels are for Kerr BHs, the bottom panels for JP BHs with ea = 10 a*. The horizontal red solid line indicates the spin 
parameter a* = 0.98 in the top left panel, and the Novikov-Thorne radiative efficiency tjnt = 0.234 (which requires a, = 0.98 
in the Kerr metric) in the right panels. From Bambi (2015c I, under the terms of the Creative Commons Attribution License. 


candidates have a mass around 10 Mq. If the compan¬ 
ion is a few Mq, the BH candidate cannot significantly 
change its mass and spin angular momentum even if it 
swallows the whole star. If the companion is heavy, its 
lifetime is very short, and it is impossible to transfer the 
necessary amount of material to spin the BH candidate 
up even in the case of accretion at the Eddington rate. In 
the end, a BH candidate cannot get more than a few Mq 
from the companion star, and for a 10 Mq object this is 
arguably not enough to appreciably change its spin pa¬ 


rameter for sufficiently large values of spin (King & Kolb 


19991. 


Moreover, while there are still uncertainties in the an¬ 
gular momentum transport mechanisms of the progeni¬ 
tors of stellar-mass BHs, it is widely accepted that the 
gravitational collapse of a massive star with solar metal- 


licity cannot create fast-rotating remnants (Woosley & 
Bloom, 2006 Yoon et al. 2006). However, this is not 


what we observe. Assuming the Kerr metric, we see 
BHs with spin close to 1, see Tab. |l] For instance, the 
BH candidate in GRS 1915-1-105 has a* > 0.98 and 


M = 12.4 ± 2.0 Mq, while the stellar companion’s mass 
is M = 0.52 ± 0.41 Mq. In the case of high-mass X-ray 
binaries, the BH candidate in Cygnus X-1 has a* > 0.98 
and M = 14.8 ± 1.0 Mq, while the stellar wind from 
the companion is not an efficient mechanism to transfer 
mass. 

A speculative possibility to explain this puzzle is to 
admit that BH candidates are not the Kerr BHs of gen¬ 


eral relativity (Bambi 2015cI. The top panels in Fig. 26 


show the evolution of the spin parameter (left panel) and 
of the Novikov-Thorne radiative efficiency (right panel) 
of a Kerr BH as a function of the accreted mass for a few 
different BH initial masses. At the beginning, the BH 
is assumed to be non-rotating. We can see that a Kerr 
BH with M ~ 12 Mq and a* > 0.98 (corresponding to 
a radiative efficiency ? 7 nt > 0.234, which is the actual 
quantity measured with the continuum-fitting method) 
had to swallow about 6 Mq. The bottom panels are as 
the top ones for a JP metric with £3 = ka^ and k = 10. 
In this case, the object needs to get a smaller mass from 
the companion, about 2 Mq, to acquire a radiative effi- 
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ciency to explain the Kerr measurement of its spin. More 
detailed calculations confirm that such an exotic scenario 


could explain this puzzle (Bambi 2015c I. If it were pos¬ 


consideration, because in some models the gravitational 
wave signal is very different from that expected in general 
relativity, while in other frameworks it is not. 


sible to estimate the mass transfer to pass from a BH 
with low radiative efficiency to a BH with high radiative 
efficiency, it could be possible to test the Kerr metric. 


VI. COMPARISON BETWEEN TESTS WITH 
ELECTROMAGNETIC AND GRAVITATIONAL WAVE 
OBSERVATIONS 

The announcement by the LIGO/Virgo Collaboration 
of the first direct detection of gravitational waves has 
opened a new window for testing gravity in the strong 


Cardoso et al. (2016) pointed out that the ringdown 


signal from a binary coalescence - like in the case of 
GW150914 - cannot be seen as conclusive proof for the 
formation of an event horizon after the merger. They 
showed that universal ringdown waveforms indicate only 
the presence of light rings, which can be possessed even 
by very compact objects without horizon. 


Konoplya & Zhidenko (2016) tried to test the Kerr 


field regime (Abbott et al. 2016a). Gravitational wave 


experiments now promise to be able to perform precise 
tests within 5-10 years. It is natural to wonder how such 
a breakthrough can affect the attempts to test BH can¬ 
didates with electromagnetic radiation discussed in this 
review paper. 

First, it is worth noting that in general it is not 
straightforward to compare the capabilities of the grav¬ 
itational wave and the electromagnetic approaches to 
test BH candidates. As already pointed out in Subsec- 
tion |H.D] they measure different things. Electromagnetic 
tests are sensitive to the motion of particles (massive par¬ 
ticles in the accretion disk and photons propagating from 
the emission point in the disk to the detection point in 
the flat faraway region). Gravitational wave tests are 
sensitive to the field equations of the gravity theory, and 
eventually one can study the evolution of perturbations 
on the background metric. 

In some alternative theories of gravity, BHs are still 
described by the Kerr metric, but the emission of gravi¬ 


metric, independently of the gravity theory, from the sig¬ 
nal of GW150914. In the coalescence of a binary, we can 
distinguish three stages: the post-Newtonian inspiral, the 
merger, and the ringdown. The post-Newtonian inspiral 
may be the same in many alternative theories of gravity 
if the metric deviates from the Kerr solution only in the 
near horizon region (this is not universally true and de¬ 
pends on the specific theory of gravity). The merger is a 
very short and complicated stage. The ringdown phase 
may be the most suitable to test strong gravity. As a fur¬ 


ther simplification, Konoplya & Zhidenko (2016) studied 


the quasi-normal frequencies of a scalar field in the de¬ 
formed background. In general relativity and in other 
theories of gravity, these frequencies are not much differ¬ 
ent from those of the gravitational waves derived from 
the field equations. However, this is not always true and 
there are also examples of gravity theories in which the 
scalar field quasi-normal frequencies can be quite differ¬ 
ent from those of the gravitational waves. The analysis 


Konoplya & Zhidenko (2016) showed that it is diffi- 


tational waves is typically different (Barausse & Sotiriou 


cult to constrain deviations from the Kerr solution from 
GW150914 because the measurement of the spin and the 
deformation parameter are correlated. This would re¬ 
main true even in the presence of the detection of addi¬ 
tional modej^ 


2008 

): in such a case, only the gravitational wave ap- 

Gardenas-Avendano et al. 

(2016b 

) and 

Bambi & Nam- 


proach can test the model. The contrary is also possible. 
The existence of non-minimal interaction terms between 
the electromagnetic and the gravitational fields or, more 
in general, the presence of new fields leading to a vio¬ 
lation of the weak equivalence principle may affect the 
motion of photons/particles without altering the emis¬ 
sion of gravitational waves. In similar frameworks, only 
the electromagnetic approach can detect deviations from 
standard predictions. 

Bearing in mind these fundamental differences be¬ 
tween the gravitational wave and the electromagnetic ap¬ 
proaches, some preliminary studies have already investi¬ 
gated the constraints that can be inferred by the LIGO 
data of GW150914 and compared with the constraining 
power of some electromagnetic techniques. 


palliwar (2016) studied the constraining power of, respec¬ 


tively, iron line and QPOs for the metric discussed in 


Konoplya & Zhidenko (2016) in order to compare the 


electromagnetic and the gravitational wave approaches. 
The conclusions of these preliminary studies can be sum¬ 
marized as follows. In the presence of high-quality data 
and the correct astrophysical model, the iron line method 
can provide strong constraints on the spin and the defor¬ 
mation parameter. The technique may thus be able to 
compete, or to be complementary to, the gravitational 
wave approach. The necessary high-quality data may be 
already available in the case of BH binaries. The system- 
atics is currently the main concern about this approach 


Yunes et al. (2016) considered a number of alterna¬ 


tive theories of gravity and showed that the data of 
GW150914 can already strongly constraint some mod¬ 
els. This depends on the specific gravity theory under 


The data of GW150914 are consistent with a single damped sinu¬ 
soid, which can be naturally interpreted as either the I = m = 2 
mode or the I = 2, m = —2 mode since the event was a binary 
BH merger. 
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and there is not a common consensus if eventually one 
can really get reliable and accurate measurements of the 
background metric. QPOs can provide very precise mea¬ 
surements, but they are not able to break the degener¬ 
acy between the spin and possible deviations from the 
Kerr metric. The parameter degeneracy cannot be bro¬ 
ken even imagining some very accurate frequency mea¬ 
surements in the future. 


VII. SUMMARY AND CONCLUSIONS 

Astrophysical BH candidates can be naturally inter¬ 
preted as the Kerr BHs of general relativity. However, a 
direct observational evidence that the spacetime geome¬ 
try around these objects is really described by the Kerr 
solution is still lacking and, at the same time, a num¬ 
ber of theoretical models suggest the possibility of new 
physics and macroscopic deviations from the Kerr back¬ 
ground. In this paper, I have reviewed current attempts 
to test the Kerr metric with electromagnetic radiation. 
The spectrum of the accretion disk, but even of stars 
orbiting close to a supermassive BH candidate, has fea¬ 
tures that can be used to study the metric around these 
compact objects and thus test the Kerr BH hypothesis. 


A. Current constraints 


Today, we do not have strong constraints on the actual 
nature of BH candidates and on the spacetime geometry 
around them. However, current observations are consis¬ 
tent with the Kerr BH hypothesis and can rule out some 
alternative scenarios. In particular: 


1. A number of observations are consistent with 
the fact that BH candidates have an event hori¬ 
zon (Broderick et al. 2009 [McClintock et al. 2004 


Narayan fc Heyir2602[ Narayan fc McClintock 


2008 Tournear et al. 2003). There is no direct 


proof, because such a proof is by the definition of 
event horizon impossible, but many (not all) alter¬ 
native scenarios in which these objects would not 
be BHs can be ruled out. 


2. Many compact objects made of exotic weakly- 
interacting matter can be ruled out because their 
spacetime cannot reproduce the characteristic low 
energy tail of the iron line profile expected in the X- 
ray spectrum of very fast-rotating BHs (]Bambi fc 


Malafarina 2013). The gravitational field around 


these exotic objects is never very strong (and in¬ 
deed they have no horizon), so photons cannot be 
strongly redshifted as is instead expected when near 
BHs. Current data show that sources with iron line 
profiles with a long low energy tail are common and 
therefore it is possible to rule out a number of ex¬ 
otic compact objects. 


3. The mass-quadrupole moment of the supermassive 
BH candidate in the quasar OJ287 has been mea¬ 
sured and it is consistent with the Kerr predic¬ 


tion at the level of 30% (Valtonen et al. 2010). 


This bound is weak, because BHs of different types 
and from different gravity theories are usually quite 
similar. However, if we consider a compact object 
like a neutron star, the value of its mass-quadrupole 
moment is a few times that of a Kerr BH with the 
same mass and spin. 

4. The data of the thermal spectrum of the accre¬ 
tion disks of the stellar-mass BH candidates in 
GRS I9I5-I-I05 and in Cygnus X-I exclude the 
possibility of large deviations from the Kerr solu¬ 


tion (Bambi, 2014a Kong et al. 2014). These ob¬ 


servations cannot constrain any kind of deforma¬ 
tion, but they can safely rule out some deforma¬ 
tions with a strong impact on the position of the 
ISCO radius. 

5. Even in the case of supermassive BH candidates 
we can exclude some large deviations with a strong 
impact on the ISCO radius. Very deformed objects 
cannot have a high radiative efficiency. At the same 
time, one can constrain the maximum value of their 
spin parameter from considerations on the spin evo¬ 
lution. The combination of these two arguments is 
a limited allowed region on the spin parameter - de¬ 
formation parameter plane: eventually, very large 
deviations from the Kerr solution are not permit¬ 


ted (Bambi 201 Ib). 


B. Prospectives for the future 

Current constraints mainly rule out the possibility that 
BH candidates are certain compact objects made of ex¬ 
otic matter and far from forming a horizon. Indeed, 
the gravitational redshift experienced by photons emitted 
close to the surface of these objects is surely very strong 
(points 1-2 in the previous subsection). Rotation does 
not make these objects as oblate as one should expect, 
for instance, in the case of a neutron star, and this is con¬ 
sistent with the fact that these objects are BHs (points 
3-5). However, these constraints are typically unable to 
distinguish BHs from different theories of gravity. It is 
definitively challenging to reach accurate and stringent 
constraints on the metric around these objects, because 
it is necessary to have a very good astrophysical model. 
The current situation of present techniques and future 
prospectives can be summarized as follows: 

I. The continuum-fitting method is probably the best 
technique available today to test BH candidates, 
but it can only be used for the stellar-mass ones. Its 
weak points are that mass, distance, and inclination 
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angle of the source must be obtained from indepen¬ 
dent measurements (with current methods, the un¬ 
certainties are usually large and systematic effects 
are possible), and that corrections to the blackbody 
spectrum depend on certain atmosphere models. 
With more precise measurements of BH masses, 
distances, and viewing angle in the future, it is pos¬ 
sible to obtain stronger constraints on the spin (if 
we assume the Kerr metric) or on the spin and pos¬ 
sible deviations from the Kerr solution. However, 
the spectrum is typically degenerate with respect 
to the spin and the deformation parameters, so it 
is impossible to test the Kerr spacetime without in¬ 
dependent measurements of the metric. This seems 
to be an intrinsic limitation of this method. 

2. The analysis of the iron Ka line is potentially quite 
a powerful tool to test the Kerr metric. It can 
be used for both stellar-mass and supermassive ob¬ 
jects. The main problem is the astrophysical model 
and the systematic effects. Currently there is no 
unanimous consensus on the possibility of using 
this technique to get accurate and reliable measure¬ 
ments of the metric around BH candidates. Some 
assumptions are still to be proven, such as the fact 
that the inner edge of the disk is at the ISCO radius 
in the hard state. The intensity prohles currently 
used are based on phenomenological models. Even 
the lamppost geometry is one among other configu¬ 
rations. In the case of supermassive BH candidates 
and with current X-ray facilities, there is an ad¬ 
ditional limitation due to the low photon number 
count in the iron line. In the presence of the cor¬ 
rect theoretical model and high quality data, the 
iron Ka line may become a leading technique to 
test general relativity. 

3. Since tests of the Kerr metric have to face the prob¬ 
lem of a strong degeneracy among the spin and 
possible deviations from the Kerr background, the 
combination of measurements sensitive to differ¬ 
ent relativistic effects seems to be the key point to 
break the parameter degeneracy and test the Kerr 
BH hypothesis. From this point of view, SgrA* may 
be one of the best sources to test the Kerr nature 
of BH candidates. A number of new observations 
may be available in the near future and the combi¬ 
nation of these measurements could test the nature 
of SgrA*. 

4. QPOs, pulsars, jet power, etc. are potential tools 
for the future. Today we cannot use these tech¬ 
niques, either because the physical mechanism is 
not well understood, or because there are no cur¬ 
rent measurements, or the measurements are not 
yet good enough. Some of these techniques may 


eventually work and be useful to test the Kerr met¬ 
ric. Other approaches may remain just as an idea. 
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